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Abstract 



First we review the definition of a negative point mass singularity. Then we examine 
the gravitational lensing effects of these singularities in isolation and with shear and 
convergence from continuous matter. We review the Inverse Mean Curvature Flow 
and use this flow to prove some new results about the mass of a singularity, the 
ADM mass of the manifold, and the capacity of the singularity. We describe some 
particular examples of these singularities that exhibit additional symmetries. 
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Chapter 1 



Introduction 



General Relativity is currently the accepted model of the physics of the universe on 
large scales. The theory principally consists of three parts. 

• The universe is modeled by a four dimensional Lorentzian manifold. 

• The geometry of this manifold is given by the Einstein equation: 

Rc> —^Rg^ = 87rT Mi ,. (1.1) 

The left hand side is called the Einstein Tensor and the right hand side is the 
stress energy tensor for the system. The stress energy tensor is a property of 
the matter in the system. 

• Particles, in the absence of other forces, move on timelike geodesies. 

The theory has had remarkable success. Its major accomplishments include explain- 
ing the orbital perihelion precession of Mercury, differing values for the bending of 
light near a massive body, and a working cosmological model running up to moments 
after the Big Bang. 
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In particular, Gravitational Lensing has provided the most useful method for 
observing and studying that portion of the matter in the universe that does not emit 
or absorb light. This is currently believed to comprise approximately 85% of the 
total matter in the universe. 

Despite its many successes, there are a number of difficulties that arise when 
working in the theory. The Einstein Equation, even in a vacuum, is a nonlinear 
second order hyperbolic differential equation in the metric. This makes exact solu- 
tions difficult to find. While the behavior of small test particles is easily given by 
the geodesic condition, the behavior of continuous masses, with features like tension, 
pressure, et cetera, are not given by the theory, but require an external derivation. 
Furthermore the relevant manifold is Lorentzian, not Riemannian, removing many 
powerful tools. 

To avoid many of these difficulties, one may study Riemannian General Rela- 
tivity. This is the study of Riemannian 3-manifolds that could arise as spacelike 
hypersurfaces in a spacetime in General Relativity. The Einstein Equation is trans- 
lated into equations about the metric on this spacelike hypersurface and its second 
fundamental form. The properties of the matter in the theory are replaced by condi- 
tions such as the Dominant and Weak Energy conditions, which are also conditions 
on the metric of the hypersurface and its second fundamental form. Many of the 
important questions in General Relativity have analogues on Riemannian General 
Relativity. For example, the Penrose Conjecture can be restricted to the Riemannian 
Penrose Inequality. 
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This thesis consists of a study of the properties of Negative Point Mass Singu- 
larities. The motivating example of which is the spacial Schwarzschild metric with 
a negative mass parameter: 



In addition to being historically and physically important, the Schwarzschild solution 
is of particular mathematical interest since it is the case of equality of the Riemannian 
Penrose conjecture, and, in the case when m — 0, it is the case of equality of the 
Riemannian Positive Mass Theorem. Thus this metric, and its generalizations, show 
promise as objects of study. 

This thesis consists of three main topics. After laying out the necessary definitions 
in Chapter [21 we examine the gravitational lensing effects of these singularities in 
Chapter [3J Next we summarize the work of Huisken and Ilmanen on Inverse Mean 
Curvature Flow in Chapter [U We then use this to prove a number of results in 
Chapter [5j In Chapter [6] we show what additional information can be gained if the 
singularities possess additional symmetries. 




m < 0. 



(1.2) 
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Chapter 2 



Definitions 



2.1 Asymptotically Flat Manifolds 

Physically, we want to make sure our manifolds represent isolated systems. A precise 
formulation of this is asymptotic flatness. 

Definition 2.1.1. Q7J) A Riemannian 3-manifold (M,g) is called asymptotically flat 

if it is the union of a compact set K, and sets Ei diffeomorphic to the complement 

of a compact set i<Q in M. 3 , where the metric of each Ei satisfies 

C C 

\g%j - < ]-p \9ij,k\ < —2 (2.1) 

\ x \ \x\ 

as |x| — > oo. Derivatives are taken in the flat metric Sij on x G M 3 . Furthermore the 



Ricci curvature must satisfy 



Rc > — |. (2.2) 

\x\ 



The set Ei is called an end of M. 



A manifold may have several ends, but most of our results will be relative to a 
single end. 
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In [T], Arnowitt, Deser and Misner define a geometric invariant that is now called 
the ADM mass. 

Definition 2.1.2. The ADM mass of an end of an asymptotically flat manifold is 

wadm = hm -— / (gij.i - gu.j) n J dfx. (2.3) 

7--+OC 107T J gr 

This quantity is finite exactly when the total scalar curvature of the chosen end is 
finite. This definition appears to be coordinate dependent, however in [1] the authors 
show that it is actually an invariant when 

/ \R\ < oo. (2.4) 

J M\K 



2.2 Quasilocal Mass Functionals 

While the ADM mass provides a definition for the total mass of a manifold, or the 
mass seen at infinity, there is no computable definition for the mass of a region. The 
two of most relevance are the (Riemannian) Hawking mass and the Bartnik mass. 

Definition 2.2.1. The Hawking mass of a surface £ is given by 



m H = 




Consider the Hawking mass of a surface, S, in R 3 , and note that H — K\ + k 2 , 
where Ki are the principal curvatures of E. Thus if 2 = k\ + n\ + 2kik, 2 - Hence, for 
a sphere in IR 3 , 

/ H 2 = / 2K + k\ + k 2 2 > / 2K + 2K = 8ttx(S) = 16vr. 
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Thus the Hawking mass of a sphere is always nonpositive in IR 3 . Furthermore the 
Hawking mass can be decreased by making the surface £ have high frequency oscil- 
lations. These two observations lead to the conclusion that the Hawking mass tends 
to underestimate the mass in a region. 

The other quasilocal mass functional of interest is the Bartnik mass defined in 

0- 

Definition 2.2.2. Let the asymptotically flat manifold (M,g) have nonnegative 
scalar curvature. Let Q be a domain in M with connected boundary. Assume M has 
no horizons (minimal spheres) outside of Q. Call any asymptotically flat manifold 
(M, g) acceptable if it has nonnegative scalar curvature, contains an isometric copy 
of Q, and has no horizons outside of Q. Then the Bartnik mass, meffi) of Q is 
defined to be the infimum of the ADM masses of these acceptable manifolds. 

The positive mass theorem guarantees that this mass will be positive if the interior 
of £ fulfills the hypotheses of the theorem. The Bartnik mass is very difficult to 
compute. The only cases where it is known are when the surface can be embedded in 
the exterior region of the Schwarzschild spacial metric or in M 3 . The Schwarzschild 
metric is the case of equality of the Riemannian Penrose inequality and M 3 is the 
case of equality for the positive mass theorem. 
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2.3 Definition and Mass of Negative Point Mass Singulari- 
ties 



The basic example of a negative point mass singularity is the negative Schwarzschild 



where m < 0. This manifold fails the requirements of the positive mass theorem since 
it is not complete: geodesies reach the sphere at r = — m/2 in finite distance. A 
straightforward calculation shows that the ADM mass of this manifold is given by m. 
Furthermore the far field deflection of geodesies is the same as for a Newtonian mass 
of m. These results are identical to the same results for a positive mass Schwarzschild 
solution. 

Two important aspects of this example will be incorporated into the definition 
of a negative point mass singularity. One is that the point itself is not included. 
To justify the use of the word "point" we have to describe the behavior of surfaces 
near the singularity. The manifold in that region should have surfaces whose areas 
converge to zero. In addition the capacity of these surfaces should go to zero. The 
second aspect is the presence of a background metric, in this case the flat metric. 
This background metric will provide a location where we can compute information 
about the singularity. 

This example motivates the following definition. 

Definition 2.3.1. Let M 3 be a smooth manifold with boundary where the boundary 
is compact. Let II be a compact connected component of the boundary of M. Let 



solution. This is the manifold M 3 \ B_ m / 2 with the metric 




(2.6) 
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the interior of M be a Riemannian manifold with smooth metric g. Suppose that, for 
any smooth family of surfaces which locally foliate a neighborhood of II, the areas 
with respect to g go to zero as the surfaces converge to II. Then II is a Negative 
Point Mass Singularity. 

A particularly useful class of these singularities are Regular Negative Point Mass 
Singularities. 

Definition 2.3.2. Let M 3 be a smooth manifold with boundary. Let the boundary of 
M consist of one compact component, II. Let (M 3 \ II, g) be a smooth Riemannian 
manifold. Suppose that, for any smooth family of surfaces which locally foliate a 
neighborhood of II, the areas with respect to g go to zero as the surfaces converge 
to II. If there is a smooth metric g on M 3 and a smooth function Tp on M with 
nonzero differential on II so that g = Tp 4 ]}, then we call II a Regular Negative Point 
Mass Singularity. We call the data (M 3 ,g, Tp) a resolution of II. 

Notice that while IT is topologically a surface, and it is a surface in the Riemannian 
manifold (M 3 ,g), the areas of surfaces near it in (M 3 \ II, g) approach zero, so we 
will sometimes speak of II as being a point p, when we are thinking in terms of the 
metric g. Furthermore, notice that the requirement that areas near II go to zero 
under g tells us that Tp — on II. 

We can define the mass of a Regular Negative Point Mass Singularity as follows: 

Definition 2.3.3. Let (M 3 ,~g,Tp) be a resolution of a regular negative point mass 
singularity p — II. Let v be the unit normal to II in g. If the capacity of p is zero, 
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then the regular mass of p is defined to be 

\ 3 / 2 



m n (p) = -~ (IJ vfflWdAjj 



(2.7) 



If the capacity of p is nonzero, then the mass of p is defined to be — oo. 

See Chapter [5] for a discussion of the capacity of points like p. We can also define 
the mass of a negative point mass singularity that may not be regular. 

Definition 2.3.4. Let (M 3 ,g) be an asymptotically flat manifold, with a negative 
point mass singularity p. Let Ej be a smooth family of surfaces converging to p. 
Define hi by 

Ahi = (2.8) 
lim hi = 1 (2.9) 

x— >oo 

hi = on E;. (2.10) 

Then the manifold (M, hfg) has a negative point mass singularity at Ej = Pi which 

is resolved by (M, g, hi). Define the mass of p to be 

I / 1 r \ 3 / 2 

sup lim — (— / u(hi)^ 3 dA ) = sup lim m-&(pi)- (2-H) 



Here the outer sup is over all possible smooth families of surfaces {E^} which converge 
to p. 

A straightforward calculation shows that if the capacity of p is non-zero, then 
the mass of p is — oo. This is the definition we will be working with. However, an 
alternative definition for the mass of a negative point mass singularity follows. 
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Definition 2.3.5. Let (M 3 ,g) be an asymptotically flat manifold, with a negative 
point mass singularity p. Choose a function h that satisfies the equations 

Ah = (2.12) 

h=- + 0[\) (2.13) 



lim/i = oo. (2.14) 

X— ¥p 

Now define surfaces E$ = {x|/i(x) = £} and functions v 5 ^) = 1 — ^A- Then the 

manifold (M, <pfg) has a negative point mass singularity at S t = p t which is resolved 

by (M, (p t ). Define the mass of p to be 

I / 1 r \ 3/2 

sup lim — I — / v((p t ) 4 ^ dA I = sup lim TOR,(pt). (2-15) 

Here the outer sup is over all possible h's. 

2.4 Fundamental Results 

Before we continue we must verify that these definitions are consistent. First it must 
be verified that the regular mass of a regular singularity is indeed intrinsic to the 
singularity, as shown in [6]. 

Lemma 2.4.1. The regular mass of a negative point mass singularity is independent 
of the resolution. 

Proof. Let (M 3 ,g, Tp) and (M 3 ,g, tp) be two resolutions of the same negative point 
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mass singularity, p. Then define A by tp = Xtp. Thus we note the following scalings: 

g = X 4 g (2.16) 

dA = X 4 dA (2.17) 
lp = \~ 1 tp (2.18) 
^=A~ 2 I7 (2.19) 

Now note that since tp, Tp — on II, II, 

= A~ 2 77 (A^) = \~*V (tp) + A- 4 77(A)^. (2.20) 

The last term, \~ A u(X)Tp, needs discussion. Both Jp and tp are smooth functions with 
zero set II and they both have nonzero differential on II. Thus A is smooth. To see 
this choose a coordinate patch on the boundary where II is given by x = 0. Then 
Taylor's formula tells us that 

'1 dip 



x _ Jo ^(xs,y,z)ds 
Jo ^{xs,y,z)ds 

which is a nonzero smooth function. Thus since tp goes to zero on II, this last term 
is zero on II. Thus the mass of p using the (M 3 ,g, tp) resolution is given by 

m R (p) = -l(±J_ ^ Vs dAj ^ (2.22) 

^ x ^ 4/3 x ™) 3/2 (2 - 23) 

=-K^ [a ~^ )]VSa4S )^ (2 ' 24) 

=-\{li v{ip)4/3M ) 3/2 ■ (2 - 25) 



□ 
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Definition 12.3.41 seems to involve the entire manifold, as the definition of hi takes 
place on the entire manifold. However that isn't the case. The mass is actually local 
to the point p. 

Lemma 2.4.2. Let (M 3 ,g) be a manifold with a negative point mass singularity p. 
Let g be a second metric on M that agrees with g in a neighborhood of p. Then the 
mass of p in (M 3 ,g) and (M 3 ,g) are equal. 

Proof. The goal is to show that for any selection of {Si}, the series rn-^ipi) and 
^r(Pi) obtained in the calculation of the mass of p, with respect to (M,g) and 
(M,g) converge. Let S be a smooth, compact, connected surface separating p from 
infinity and contained in the region where g and g agree. Fix i large enough so that 
Ej is inside of S, and suppress the index i on all our functions. Then define the 
functions h, h by 

h = h = on Sj 
lim h = lim h = 1 

Ah = Ah = 0. 

Here A and A denote the Laplacian with respect to g and g respectively. 

Now inside S, A = A since g = g. Thus there is only one notion of harmonic, 
and h and h differ only by their boundary values on S. Let e = 1 — min^j/i, h}. 
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Consider the following two functions / and / + defined between S and 

/" = /+ = on E, 

Af- = A/+ = 

f~ = 1 — e on S 

f + = 1 on 5. 

Thus by the maximum principle, we have the following inside S 

f + >h,h>r. (2.26) 

Furthermore, since all four functions are zero on Ej, 

> v(h),u(h) > u(f-). (2.27) 

Here v is the normal derivative on Ej. Now define £(<p) by the formula 

£{<p) = [ v(<p) i/3 dA. (2.28) 

Then the ordering of the derivatives gives the ordering 

£(f + )>£(h),£(h)>£(f-). (2.29) 

However, since f~ = (1 — e)/ + , 

!/(/-) = (l-c)i/(/ + ), (2.30) 

hence, 

£(/") = (l-e) 4 ^(/+). (2.31) 

Now, without loss of generality assume that the limit of the capacities of {Ej} is 
zero, as the mass would be — oo otherwise. Thus as i — > oo, Ej has capacity going 
to zero. Hence ej goes to zero, and so £(f^)/£(f^) goes to 1. Thus equation (12. 29[) 
forces £{hj) and £{hj) to equality. This forces the masses of pt in the two metrics to 
equality as well. □ 
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Corollary 2.4.3. In Definition \2.3.4\ we may replace the condition that (fi be one at 



infinity with the condition that (fi be one on a fixed surface outside Ej fori sufficiently 
large. 

This mass also agrees with the regular mass when the singularity is regular. 

Lemma 2.4.4 ([3J). Let (M 3 ,g) be an asymptotically flat manifold with negative 
point mass singularity p. Let p have a resolution (M 3 ,g,ip). Then the regular mass 
of p equals the general mass of p. 
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Chapter 3 



Gravitational Lensing by Negative Point 
Mass Singularities 



3.1 Gravitational Lensing Background 

One of the first testable predictions of general relativity was the difference in the 
deflection of light by gravity. This effect was first confirmed during the 1919 solar 
eclipse. Since then gravitational lensing has become an powerful tool for astronomy 
in general and cosmology in particular. Gravitational lensing has made it possible 
to detect the presence of dark matter by observing its effects on background images. 

In this chapter we will develop the properties of lensing by negative point mass 
singularities in the setting of accepted cosmology. We will make a number of as- 
sumptions based on that cosmology that will allow us to obtain simple formulas for 
gravitational lensing by negative point mass singularities. Then we will character- 
ize their lensing effects and compare them to positive mass point sources. We will 
not cover the entire field of gravitational lensing, but only develop enough for our 
purposes. 

15 



We will restrict ourselves to negative point mass singularities that agree with a 
negative Schwarzschild solution to first order. We will find that the lensing effects of 
these singularities in the presence of continuous matter and shear can be duplicated 
by configurations with positive mass lenses. 

We will follow the presentation given in [9] . We will differ from this presentation 
by using geometrized units where c — G — 1. We will also consider lens potentials 
outside of the scope of [9]. 

3.2 Cosmology for Gravitational Lensing 

To simplify the calculations involved in lensing, we will make use of a number of 
assumptions about the configuration of our system and its behavior. These assump- 
tions are based on the scales and phenomenology of astronomy. Our first assumption 
is one of cosmology. 

Assumption 3.2.1. The universe is described by a Friedmann-Robertson- Walker 
cosmology. 

The Friedmann-Robertson- Walker model is an isotropic homogeneous cosmology 
filled with perfect dust. We will not develop this cosmology from these properties 
but merely take it as a given. In this cosmology, the universe is modeled as a warped 
product with leaves K and fibers given by either M 3 , H 3 or S 3 . Thus we have the 
metric 

ds 2 = -dr 2 + a 2 (r)dS 2 K . (3.1) 

Where Sk is H 3 , R 3 , or S 3 when K = —1,0,1, respectively. The coordinate r is 
time as measured by the isotropic observers. This is called "cosmological" time. The 
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function a(r) gives the scale of the universe, and has the following relationship to r 
depending on K 

A A 

K — l a = — (1 — cos-u) t = — (1 — sin it) (3.2) 

K = a=( 9 4y 3 r 2/3 (3.3) 



A A 
K = — 1 a = — (cosh-u — 1) t = — (sinhit — u) . (3.4) 

2 2 



The metric on the fibers is given by 

dR 2 
1-KR 2 

If we write 



dS 2 K = dK + R 2 (d9 2 + sin 2 9 d V 2 ) . (3.5) 



sin x if X = 1 

n = »m l< (\) = { X if* = (3.6) 
sinhx ifx = -l, 

then we can rewrite the fiber metric as 

dS 2 K = d X 2 + sin^ X {dO 2 + sin 2 6 d V 2 ) . (3.7) 

The differences between all of these metrics are small on scales small compared to 

the size of the universe, where most of our calculations will take place. We also 
introduce a second time coordinate, t, by the following equation 

(3-8) 

We can use this to rewrite the metric on the universe as 

ds 2 (t) =a 2 (t) (-dt 2 + dS 2 K ) . (3.9) 

For this reason t is called "conformal time." 
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In our discussion of the geometry of a lens system, we will need to discuss the 
distances between the observer, lens and source. There are a number of options 
available, depending on which equations one wants to simplify. We will use "angular 
diameter distance." This distance is defined by the ratio the physical size of an 
object to the angular size of the object seen by the observer. We denote the distance 
from the observe to the lens by dz,, the distance from lens to source by d^s, and the 
distance from observer to source by ds- For lensing on small scales (nearby galaxies) 
the universe is almost flat, and since the bending angle is generally small, we can use 
d s ^ d LtS + d L . 

We will also need the redshift. As the universe expands, the wavelength of light 

from distant sources is lengthened. This is equivalent to clocks appearing to running 

more slowly. This is properly associated to an event, but since it changes slowly 

compared to the size and duration of a lensing event, we will be assuming it is 

locally constant. We can define the redshift of a time, r, (or an event) by 

*=^-l. (3.10) 
a(r) 

3.3 Geometry of Lens System 

We make a number of assumptions about the geometry of a gravitational lensing 
system. First we assume that the lens isn't changing on the time scale it takes for 
light to cross the lens. This is valid since most objects evolve at speeds much less 
then the speed of light. We can also assume that the lens is stationary. Any relative 
motion will be attributed to the source. 

Assumption 3.3.1. The geometry of the spacetime is assumed to be unchanging 
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on the time scale of the lensing event. 

Furthermore, since almost all sources are "weak" we also assume that the en- 
tire system lies in the weak regime, where can approximate the metric by a time- 
independent Newtonian potential, ip, given by 

(p(x) = -a 2 L f P(X l dx. (3.11) 
Jr 3 \\ x ~ x II 

Here cll is the value of a when the light ray is interacting with the lens. 

Assumption 3.3.2. The spacetime metric of gravitational lens system is given by 
ds 2 = -(1 + 2<p) dr 2 + a 2 (r)(l - 2(p) dS 2 K (3.12) 

= a 2 (t) [-(1 + 2y?) dt 2 + (1 - 2ip) dS 2 K ] . (3.13) 

Furthermore, tp, the Newtonian potential, is much smaller then unity. We will 
also assume that while the light ray is interacting with the lens, a(t) is constant. 
Furthermore, since on the scale of this interaction, the universe is almost flat, we 
will assume K = during the interaction. Thus near the lens we have the following. 

Assumption 3.3.3. During the interaction of the light ray with the lens, the metric 

of the lens can be assumed to be 

ds 2 L = a\ [-(1 + 2y?) dt 2 + (1 - 2<p) dS 2 ] . (3.14) 

Furthermore, dt can be approximated by -^-dr. 

We choose dimensionless Euclidean coordinates, (x 1 ,x 2 ,x 3 ), centered at the lens 
so that dS 2 = 5ij and so that the line of sight to the lens is along the x% axis. We will 
also use proper coordinates (ri,r 2 ,C) = a-L • (^1^2,^3), and will use the coordinates 
r = (ri,r 2 ) on the lens plane. We will use proper coordinates s = (s 1 ,s 2 ) in the 
source plane. 
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3.3.1 Mass Densities, Bending Angle and Index of Refraction 

Since we are in the static weak field limit, the Einstein equation reduces to the time 
independent Poisson equation for ip. Thus we get the three dimensional Poisson 
equation 

A 3D <p(x) = Aira\p(x). (3.15) 

Where p is the density above the background of the lens. The solution to this 
equation is 

" w = - ai f«£^r (316) 

JR \\ x x II 

Since we are assuming that the lens is planar, it is useful to project the three dimen- 
sional potential and density into the lens plane. Integrating p along the line of sight 
gives us the surface mass density of the lens cr(r). This integral is really only from 
—dL,s to di, but since we are assuming that p is zero except near the lens we can 
extend this integral to (—00, 00). Thus we get that 

a(r)= [ pir^QdC. (3.17) 



Integrating the three dimensional Poisson equation along the ( axis gives us 

= gf + g + g)« = A»^, «. (3,8) 

The term integrates to zero since tp is zero when ( = ±00. If we define the 
surface potential of the lens by 

tt(r) = 2 / p(ri,r 2 ,C)dC, (3-19) 
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then \1/ satisfies the two dimensional Poisson equation 

A 2D ^(r) = 8w(r). (3.20) 



This is solved by 



ij/(r) = 4 / a(r') In 
Jm 2 



dr', (3.21) 



d 

for any constant d . We will generally choose d = d L . 

It will become useful to think of the potential of the lens giving a refractive index 

to the spacetime. The index of refraction of a medium is the reciprocal of the velocity 

of light in that medium. We want to calculate the velocity of light in our metric, 

ds 2 L , relative to the flat metric given by a 2 L 5ij. The velocity of light in the metric 

ds 2 = -A{x)dt 2 + B{x)dS 2 (3.22) 



is given by the ratio 



n=\h- ( 3 - 23 ) 



In the metric ds 2 L we get 



" = ^irf (3 ' 24) 



to first order in (p. Since p> is much smaller then unity, we can ignore higher order 
terms. 

Now we can look at the bending angle of the lens. We approximate the light ray 

from the source to the observer by a broken null geodesic with the break at the lens 

plane. We compress all the bending in the ray due to the lens into this corner of the 

light ray. Call the tangent to the incoming light ray Tj(r) and tangent to the final 

ray Tf(r). Then we define the bending angle by 

a{r)=T f {r)-T l {r). (3.25) 
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Here we have parametrized the vectors by the impact parameter r. Continuing 

with the standard geometric optics approximation, we parametrize the spatial path 

R(s) = (Rx(s), R 2 (s), R 3 (s)) of the light ray by arclength, s, in the background 

metric (ds 2 = a 2 L Sij). Then light rays are characterized by the equation 

d ( dR\ „ 

' n — = Vn. (3.26) 



ds \ ds J 

Here V is the flat gradient. Now define the quantities T = ^ and K = ^ as the 
tangent and curvature vectors of the curve R(s). Plugging these into equation f !3.26|) 
gives the equation 

{Tn)T + nK = Vn. (3.27) 

Since K is perpendicular to T the transverse gradient is given by 

V ± n = nK. (3.28) 

Solving this for K gives us 

K = ~ (-2V±(p) (1 - 2v?) _1 ~ -2V ± y? (3.29) 
n 

to first order in (p. Since this angle is small, the light rays are almost perpendicular 
to the lens plane, so we can replace Vj_ with the gradient in the (n,^) plane, V r , 
and we can integrate over ( to find the total K for the entire light ray. This total K 
tells us how far the tangent vector has turned, hence the bending angle is 

a(r) = 2 J VMn,r 2 ,0dC (3.30) 

Pushing the integral inside the gradient gives us 

a(r) = W(r). (3.31) 
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3.3.2 Fermat's Principle and Time Delays 

One could use equation (13.311) to try and work out the effect of a lens, but instead 
we will follow [9 J and use the following principle. 

Proposition 3.3.4 (Fermat's Principle). A light ray from a source (an event) to an 
observer (a timelike curve) follows a path that is a stationary value of the arrival 
time functional, T , on paths. 

Here we are only considering paths v r , that are broken geodesies from the source 
S, to the observer O, parametrized by the impact parameter r = (r 1; r 2 ) where the 
ray crosses the lens plane. For a given source location s in the source plane, we look 
at the time delay function T s (r) that gives the time delay for a light ray that goes 
from s to r, bends at r, and then continues on to O. Technically the time delay 
function and the arrival time function differ by some reference value. That reference 
value is the time the ray would have taken in the absence of the lens. We denote 
this unlensed references path by uq. Fermat's principle tells us that the images of a 
source at s are given by solutions to the equation 



To calculate this we will separate the time delay into two components. One is the 
effect of the longer path the light ray takes. This is called the geometric time delay, 
T g . We will drop the s for the moment. The other effect is due to time passing more 
slowly in a gravitational potential as seen by a distant observer. This is called the 
potential time delay, T p . The travel time for the unlensed ray is given by 



V r T s (r) = 0. 



(3.32) 




(3.33) 
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Where dl is given by the metric dS^. The travel time for the lensed ray is given by 

a L n L dl. (3.34) 



Hence the time delay is given by 

T L (r) = I ain L dl — I a L dl. (3.35) 



"0 



Here the L attached to T denotes the fact that these time delays are being measured 
at the lens plane. We will have to account for the redshift, Zl, of the lens. Thus we 
define the geometric and potential time delays by 

T p L = I a L (n L - l)dl T g L = a L (j dl - j dl) . (3.36) 

The potential time delay is easiest to calculate. Since the bending angle is small, we 
can approximate v r by u . Using equations (I3.24p and (I3.19P we can calculate the 
potential time delay as 

T/(r) = -¥(r). (3.37) 

The geometric time delay is more complicated. We will first calculate the geometric 
time delay assuming that K — 0, since we are calculating many quantities to first 
order, the geometric time delay will the same for K — ±1. For a detailed treatment 
of K = ±1, see [9]. 

First we define the dimensionless lengths Is, II, and 1^$ as the lengths of the 
spatial projections of u , and the parts of v r between the lens and observer and 
observer and source respectively. These lengths are measured in the metric dS^. 
Thus the geometric time delay measured at the lens is given by 

T g L = a L {l L + l L)S - l s ) . (3.38) 
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The law of cosines tells us 



l S = l L + l LS- 2l Lk,S COs(7T - Of) 



(3.39) 



We can approximate cos(7r — at) by —1 + a 2 /2 since a is small. This gives us 

l 2 s-ll + lis + 2Wl,s(! - «V2) (3.40) 



= (Jl + h,sf ~ hh,s® 2 - 

Isolating the term with a and factoring gives us 

(Jl + Il,s ~ Is) (h + k,s + Is) - hhji& 

IlIl.s 



' g ~ 



h + h.s + I 



-a 2 ~ 



IlIl,s - 2 
— -r—oc. 
2lo 



(3.41) 

(3.42) 
(3.43) 



We can replace II and Il,s by cllcIl, and aidL,s respectively. We would like to remove 
the reference to at. To do that we first construct the point s' in the source plane. It 
is the location that would produce an image at r in the absence of the lens. Then 
using similar triangles we compute 



a d 



L,S 



r s 
d L d s 



ds. 



(3.44) 



Plugging this into equation (13.431) gives us 

1 d L d s 



cll 2d 



L.S 



r 

d L 



s 

d s 



(3.45) 



Adding the two parts of the time delay together, and correcting for the redshift by 
multiplying by 1 + Z£ gives us 



%{r) = {l + z L ) 



d L d s 


i 


r s 


d L ,s 




d L d s 



d 



L,S 



d L d s 



(3.46) 



Taking the gradient by r and solving for s gives the lens equation 

s = —r - d LjS a{r). 
d L 



(3.47) 



25 



Viewing s as a function of r results in the lensing map. This map goes from the 

image plane to the source plane. It answers the question: "Where would a source 

have to be create an image at this location?" To further simplify this equation, we 

will nondimensionalize by introducing the following variables 

x = r/d L y = s/d s (3.48) 

^( x ) = f^f] #(r) a(x) = %£a(r) (3.49) 
\d L d s J d s 

k(x) = where cr c = f^- — . (3.50) 

a c 2ixd L d L>s 

With these variables, the lensing map becomes 

y = rj(x) = x — a(x) (3.51) 

with a = V^>. 

3.3.3 Magnification 

In addition to changing the location of the image of a source, gravitational lensing 
can also change the apparent size of an object. Since all sources are not truly point 
sources, we can really consider how a region, R, around a point s in the source plane 
gets deformed. In particular, the signed area of the image of R will be determined by 
the integral of the determinant of di]^ 1 . Due to the Brightness Theorem the apparent 
surface brightness of an object is invariant for all observers. For instance if one were 
twice as far from the sun, the total light received would be reduced by four, as would 
the area of the sun. So the observed surface brightness would remain constant. Thus 
the total light received at the observer from an extended source is scaled the same 
as the areas. See [11] for more information. 
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However, while many sources aren't point sources they are point-like. Thus the 
only effect of the magnification is to increase (or decrease) the brightness of the 
source. Thus the magnification of a point source at the location y is given by 



Locations, x, in the lens plane where this is infinite are called critical points. The 
corresponding locations r)(x) = y in source plane are called caustics. Sources on 
different sides of caustics typically have two fewer or less images then each other. It 
is often useful to consider a source that is moving in the source plane. As this source 
crosses the caustic, two of its images will increase in brightness and merge, then 
disappear, or the reverse depending on the direction in which the source crosses the 
caustic. The changing magnification of a source as it moves in the source plane is a 
useful observable called the light curve. For these curves, we add up the magnification 
of all the images, since sometimes the images are too close to be resolved. 

3.4 Lensing Map for Isolated Negative Point Mass Singular- 



The general framework we have established for gravitational lensing requires a po- 
tential function to plug into the metric in assumption 13.3.31 and following formulas. 
We will be studying singularities that agree to first order with negative Schwarzschild 
solutions. This is summarized in the following assumption: 

Assumption 3.4.1. The two dimensional surface mass density for a negative point 
mass singularity is given by 




(3.52) 



ities 




(3.53) 
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To first order, this assumption gives the correct behavior of a negative mass 
Schwarzschild solution in the weak field regime. We will restrict ourselves to cases 
that agree with this case to first order^. 

We nondimensionalize our potential by defining 

M 

m = — =—. (3.54) 
ixd\p c 

Which gives us the dimensionless surface mass density 

k(x) = 7fm5(x). (3.55) 

This gives us the dimensionless surface potential 

V>(z) =mln(||x||), (3.56) 

and the dimensionless bending angle 

a(x) = m X 2 . (3.57) 
||x|| 

Thus the dimensionless lensing map is given by 

1 ~ WJ ' (3 ' 58) 

In this case, we can solve this equation exactly to find the images of a source at y. 



x± = 7: ( \\v\\ ± VIMI 2 + 4m ) v- ( 3 - 59 ) 



Here y is the unit vector in the direction of y from the origin. As long as ||y|| > 
we get two images for each source. Both images are on the same side of the 



1 For a more detailed assumption of the positive mass analogue of this assumption see the discussion 
of point masses on p. 101 in [9]. 
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lens. If we had m > 0, then our images would be on opposite sides of the lens. The 
derivative of the lensing map is given by 

2xiX2m 




j m _|_ 2x|m 



(3.60) 



Hence the magnification is given by 



1 Itall 4 
= — = „ 4 „ • ( 3 - 61 ) 



ixir — m 2 



For X- this number is negative,so to get the total magnification we take the difference 
between signed magnifications^ 

II y II 2 — ///" 

Aitot(y) = - My = j (3-62) 

Ibll ylbll 2 + 4m 

3.5 Light Curves for Isolated Negative Point Mass Singular- 
ities 

As we calculated, the lensing map is given by 

y = v (x) = x ( 1 - J • (3.63) 



This has inverse 



x± = ^(\\y\\±J\\y\\ 2 +4m)y. (3.64) 



The inverse map tells us where an image will appear for a source located at y in the 
source plane. Since the radical is imaginary for \\y\\ < 2sJ— m, these sources aren't 
visible at all. For sources outside this disk, we get two images: x±. The reversed 



2 See Appendix I A. 1 1 for the calculation. 
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image x_ is closer to the center of the lens. For large \\y\\, the x_ image is closer and 

closer to the center of the lens, while the x + image is closer and closer to its unlensed 
location. As y gets closer and closer to the caustic \\y\\ = 2a/— m, the two images 
come together at x — y/—m. Looking at the individual magnifications, we see that 
for large \\y\\, the positive image has magnification about 1, while the negative image 
has magnification about —m 2 / \\y\\ A - As y gets closer and closer to the caustic, the 
magnification of each image increases, and is formally infinite when y = 2yJ—m and 
x± = y/—m. This curve is the critical curve. When y is inside the caustic, it creates 
no image. 

In cases where the source, lens and observer are moving relative to each other, the 
total magnification changes. We look at the total magnification, since the individual 
images are unresolvable, and the only effect of the lensing is the magnification. By 
the symmetry of the source, these paths are characterized by impact parameter, the 
distance of closest approach to y = 0. If this distance is d, and we assume that the 
source is moving at constant unit speed, the total magnification as a function of time 
is 



Figure 13.11 shows several possible paths for a moving source. The sources in these 
curves have impact parameters varying from zero to twice 2\J—m. Figure [3721 shows 
the corresponding light curves. The light curves for the sources passing inside of 
2y/—m are distinctive, but those for the source passing outside or just along the 
caustic are not. The light curve formula, equation (I3.65jl . is the same as that for a 




d 2 + t 2 + 2m 



(3.65) 



Vd 2 + t 2 Vd 2 + t 2 + 4m' 
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positive point mass. Thus if we define 



m = — m d = Vd 2 + 4m, (3.66) 

then the light curve for a source passing within d of the line of sight of a mass of m 
is the same as that for a source passing within d of the line of sight of a mass of m 

U\ d 2 + t 2 + 2m 

/*m,d(*) = n* ( 3 - 67 ) 

Vd 2 + tWd 2 + t 2 + 4m 



d 2 + 4m + t 2 - 2m 



(3.68) 



\/d 2 + 4m + tVcP + 4m + t 2 - 4m 
= /w(t). (3-69) 

Furthermore, we will see that if we introduce continuous matter we can reproduce 
the entire lensing map with a positive mass singularity. 

3.6 Complex Formulation 

Before incorporating additional features into our lens, it is helpful to reframe the 
structure of the lensing map as a map C — > C, rather than a map M 2 — >■ M 2 . First 
we will consider x as a complex number x\ + 1x2, and likewise y and 77. The lens 
equation becomes 

rj = 771 + irj 2 = [xi- 77— )+i\x 2 - 77— (3.70) 



dxi ) \ dx 

Taking complex derivatives of 77 we get 



drt \ ( drt drt 
— %■ 



dz 2 \dx\ dx- 



1 ( drji .dr] 2 .dvi . drj 2 \ 1 / drji drj 2 } _ 

+ l i^- l izr + izr) = «\-z- + -z-), (3-71) 



2 9xi dx 2 J 2 (9x2 
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which is real. Here we used the equality of the mixed partials of ip to cancel the 

imaginary terms. Differentiating with respect to ~z gives 

f / drji dr)2 



dri 1 / dri dri 

+ 1 



dz 2 \dx\ dx^ 



2 V dx\ dx 2 



+ ^ 



dr] 2 chn 
dx\ 3x2 



(3.72) 



Here no such cancellation occurs. We can also rewrite J = det {drj) as 
drji dr] 2 drji dr) 2 



J 



dx\ dx 2 0x2 dx\ 
1 drji 2 1 drji dr] 2 1 dr] 2 ' 



+ 



+ 



4 dxi 2 dxi dx 2 4 dx 



1 dr/i 2 1 drji dr] 2 1 dr] 2 ' 
4 dxi 2 dx\ dx 2 4 dx 2 



1 dr/2 2 1 dr/2 drji 1 drji 2 



+ 



+ 



4 <9xi 2 <9xi <9x 2 4 <9x 2 



drj 


2 


<9?? 


dz 




<9l 



(3.73) 



Here we extensively used the fact that I 22 - = f^. Our critical points are located 



dx 2 



where J = 0. As we noted above |^ is real so the solutions to J = look like 



dr] 

dz 



dz 



up 



(3.74) 



for some angle <p. Thus our critical curves will be curves parametrized by (p. 

The caustics will be the images of these critical curves under rj. Any points where 
the caustics aren't smooth are characterized by 

J(x) = V z (v) = 0. (3.75) 



Here Z 



dJ , • 8J 

8x2 dxi 



2i 



dJ 



and 



Z— Z— 

dz dz 



(3.76) 



To find the cusps on the caustics we just find the appropriate phase (p to solve (13. 75 p . 
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3.7 Lensing by Negative Point Mass Singularities with Con- 
tinuous Matter and Shear 



Most lensing events on the scale of individual stars take place within a host galaxy. 
In these cases, the star itself isn't the only source of distortion. Two other non-local 
factors are also important. 

Continuous matter is the first. The presence of evenly dispersed matter in the area 
of the lens can produce convergence or divergence. This enters into the dimensionless 
potential via a term like 

ipcm{x) = ^ INI 2 . (3.77) 

This gives us a lensing map of 

7)(x) = y = (1 — k)x. (3.78) 

This is clearly compatible with the complex formulation. The dimensionless mass 
density of k corresponds to a surface mass density of <j c ixk. 

The other factor is shear from infinity. The presence of a large mass nearby, such 
as a nearby galaxy, or an asymmetric distribution such as the disk of the host galaxy, 
can introduce a potential of the form 

^ sh (a;) = -- \{x\ -x£)cos20 + 2a;ia; 2 sin20] . (3.79) 

The parameter 7 defines the magnitude of the shear, and 20 determines the preferred 
direction of the asymmetric mass distribution or the direction to the large mass. By 
the symmetry of our lens, we can assume that = 0. This gives us a lensing map of 

y=\ l V 1 ° J*- (3-80) 
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In complex form this is 



y = z + 72. 



(3.81) 



Now we will incorporate all of these features into a single lens. The lensing map 
of a negative point mass lens is given by 



y = x 1 



x 



|2 ) - 



(3.82) 



In complex form this is 



y = z- 



m 



(3.83) 



Adding these potentials to that of our singularity gives us the combined lens equation 



T](x) 



1 - K + 7 

1 - k-7 



m 

x 



2""' 



(3.84) 



or 



77 = (1 — k)z + •yz — —. 

2- 



(3.85) 



Our Jacobian is 



m\ 2 



:i-k) 



(3.86) 



So to look for critical points we use equation (I3.74p with our r\ to give us 

7 + 5 = 11-*!**. 



(3.87) 



If k = 1, then we are looking for points that solve 

m 



z- 



7- 



(3.88) 



Which are just two points along the x axis. So we can assume that k 7^ 1. It 



simplifies calculation to remove k from the calculation by defining 



7* 



7 

II - k\ 



rn 
II - k\ 



e K = sgnfl - k). 



(3.89) 
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With those substitutions we now solve 




(3.90) 



This has solutions 



z±(<p) = ± 




(3.91) 



At this point the only effect that the negative sign on m has had is to rotate our 
critical curves by ir/2. 



Note that these curves are independent of e K . For 7* <g 1, the critical curve is given 
by an oval with long axis in the x^ direction. As 7* grows to 1, the critical curve 
develops a waist. For 7* > 1, we have two critical curves, small loops on the x\ 
axis. As 7* grows they shrink to the two points for k — 1. For 7* very close to 
1, both varieties of critical curves grow to infinity. When 7* = 1 our critical curves 
degenerate into two curves asymptotic to the lines X2 = ±x\. As 7* passes 1, the ends 
of the curve open up, pass through infinity, and reseal re-paired. See Figure I3T31 for 
the shapes of the critical curves for various values of shear and convergence. These 
curves are the same as in the positive mass case, rotated a quarter turn. 

To find the cusps we plug our lensing map into (13. 75 p . In our case we have 




(3.92) 




(3.93) 



and equation (13.751) is 




(3.94) 
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Figure 3.3: Critical Curves for Negative Point Mass with Continuous Matter and 
Shear 

These pictures are shown in order of increasing k, with constant 7 = 0.2. The first 
four correspond to e K — 1 and increasing 7*. At k = 1 we have two points as the 
critical set. The last four pictures correspond to e K = — 1 and decreasing 7*. Since 
the curves are independent of e K , they are the same when one replaces k by 2 — n. 
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Shear 


e« = 1 

-^cusps 


e K = -1 

Vi -^cusps 


< 7* < 3/4 





<?1, V?2 4 


3/4 < 7* < 1 


<P3, V?4, <^5, V?6 8 


¥>i, y?2 4 


Kit 


y?i, <^ 4 , y?e 6 


V?2, V?3, V?5 6 



Table 3.1: Numbers of cusps on caustics for various shear and continuous matter 
values. 



Removing the 1 — k as before, we get 



= Z; + r,.. ( 7* + Z*. 



(3.95) 



Where Z* has the same formula as Z replacing 7 and m with 7* and m*. In Ap 
pendix IA.2I we find the possible roots of this equation to be 

'3± v^r^ 



¥»i = 0, y? 2 = 7T, ^3,4 = acos 



47* 



y? 5 = 27T - (p 3l (p 6 = 27T - <p 4 . 

(3.96) 



For a given value of 7* and m*, each of these is only a solution for either e K = 1 or 
e K = — 1. The numbers of cusps for various values of 7* and e K are given in Table I3TT1 
When the number of cusps changes we expect to have higher order caustics. Looking 
at Figure E31 we see that we have four swallow tails between panels (a) and (b), 
we also have two elliptic umbilics between panels (d) and (/). The changes in the 
caustic structure between panels (b) & (d) and (/) & (h) occur at infinity. 

In terms of the number of images, we always get two images outside of the caustics. 
Inside the caustic in panel (a) we have no images, just like in the case without shear 
or convergence. In panel (b) we have four images inside the swallow tails, and none 
inside the rectangular region. In the rest we have four images inside the caustics. 

Both the critical curves and the caustics are the same as what one would get 
if one took the positive mass case, and switched k with 2 — k, and rotated all the 
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(a) 



(g) 



k = 0.75 







(d) 


k = 0.85 




k = 1.20 



X 



(b) 



(e) 



00 



k = 0.79 



k = 0.99 




k = 1.23 



(c) 



(f) 



(i) 



k = 0.80 



K = 1.15 




K = 1.50 



Figure 3.4: Caustics for Negative Point Mass with Continuous Matter and Shear 

These pictures are shown in order of increasing k with constant 7 = 0.2. The first 
four correspond to e K — 1 and increasing 7*. At k = 1 we have two points as the 
caustic set. The last four pictures correspond to e K = — 1 and decreasing 7*. 
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images by a quarter turn. 



Chapter 4 



Inverse Mean Curvature Flow 



This chapter lays out what we need of the weak inverse mean curvature flow as de- 
veloped in [7]. We will follow their exposition closely, omitting many of the technical 
details. 

4.1 Classical Formulation 

Let N be a the smooth boundary of a region in the smooth Riemannian manifold 
M. A classical solution of the inverse mean curvature flow is a smooth family x : 
N x [0, T] — > M of hypersufaces N t = x(N, t) satisfying the evolution equation 

% = -■ (4-1) 
dt H K J 

Here v is the outward pointing normal to N t and H is the mean curvature of N, 
which must be positive. We use this flow to explore the geometry near a singularity. 
For instance we will use the fact that under this flow the Hawking mass is non- 
decreasing. This result can be readily shown if we assume that the flow doesn't have 
any discontinuities or singularities. However, easy counterexamples illustrate that 
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this is overly optimistic. The simplest counterexample is given by a thin torus in 
M 3 . Such a torus has positive mean curvature approximately that of a cylinder of 
the small radius. The inverse mean curvature flow will tend to increase the small 
radius of this torus. However, if it continued without singularities or jumps, the hole 
in the torus would eventually shrink to the point where the area of the torus inside 
the hole has zero mean curvature and hence the flow couldn't be continued. Thus 
the classical flow is insufficient for our needs. 

To remedy this problem, we will follow [7] and recast the flow first in a level 
set formulation and then we will move to a weak solution. This will allow the flow 
to jump to avoid situations where the curvature of the surface would drop to zero. 
In the previous example, the flow would close the interior of the torus as soon as 
it is favorable in terms of a certain energy functional. Even with these jumps the 
Hawking mass of our surface is still non-decreasing. 

4.2 Weak Formulation 

First we establish some notation. Let (M, g) be the ambient manifold. Let h be the 
induced metric on N. Let A^j = (V^v, ej) be the second fundamental form of N. 
Then H is the trace of A with h, and H = —vH is the mean curvature vector. Let 
E be the open region bounded by N. 

The first step toward the weak formulation is a level set formulation. We assume 
that the flow is given by the level sets of a function u : M — > K. This u is related to 
our previous data by 

Ek := {x | u(x) <t}, N t := dE t . (4.2) 
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We will also need the following sets 

E+ := int {x \u(x)>t}, N+ := 0E+ (4.3) 

Where Vu ^ 0, E t = Ef and N t = N t + . 

Anywhere that u is smooth and Vu 7^ 0, then we have a foliation by smooth 
surfaces Nt, with normal vector v = Vtt/ \Vu\. The mean curvature of these surfaces 
is given by divTv(^) and the flow velocity is given by l/|Vu|, so equation (14.1 j) 
becomes 

div M C^^j = |Vw| • (4.4) 
This equation is degenerate elliptic. To remedy this we introduce the functional 

Ju{v) = jK{v)= [ \Vv\ +v\Vu\. (4.5) 
Jk 

Where K is a compact set in M. If we take the Euler-Lagrange equation of this 
functional, and replace v with u, we get back equation (14.41) . For each u we get a 
different J u . Hence what we want is a function u which minimizes its own J u . 

Definition 4.2.1. Let u be a locally Lipschitz function on the open set Q in M. 
Then u is a weak ( sub-, super- ) solution of equation (14.41) on Q exactly when 

Ju {u) < J?(v) (4.6) 

for all locally Lipschitz functions v (< u, > u) which only differ from u inside a 
compact set K contained in Q. 

It is worth noting that J„(min(t>, w)) + J u (max(v,w)) = J u (v) + J u {w). To see 
this, construct K v — {x G K \ v(x) < w(x)}, and divide the integrals on the left into 
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integrals over K v and K \ K v . Then regroup then and recombine to get the right 
hand side. This tells us that if u is a both a weak supersolution and subsolution, 
then it is a weak solution. 

We will also need a related functional of the level sets. 

Definition 4.2.2. If F is a set of locally finite perimeter, and d*F is its reduced 
boundary, then we define 



J FnK 

For any locally Lipschitz function u and compact K contained in A. We say that E 
minimizes J u in A (on the inside, outside) if 



for all F that differs from E in some compact K in A (with E C F, E D F.) A 
similar argument tells us that if E minimizes J u exactly when it minimizes J u on the 
inside and outside. 

These two formulations are equivalent. 

Lemma 4.2.3. Let u be a locally Lipschitz function in the open set Q. Then u is a 
weak (sub-, super-) solution of equation (I4.4p exactly when for each t, E t = {u < t} 
minimizes J u in Q (on the inside, outside). 

Proof. Lemma 1.1 in [7]. □ 

We now define the initial value problem. Let Eq be an open set with C 1 boundary. 
We say that u G Cj°' c and the associated E t for t > is a weak solution of f)4.4p with 




(4.7) 



J«(E) < J U (F) 



(4.8) 
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initial condition E if either 



Eq = {u < 0} and u minimizes J u on M \ E Q 



or 



(4.9) 



E t = {u < t} minimizes J u in M \ Eq for t > 0. 



These two conditions are equivalent by Lemma 1.2 in [7]. Showing the regularity 
of Nt and is nontrivial, but we won't reproduce it here. 

Theorem 4.2.4. Let n < 8. Let U be an open set in a domain Q. Let f be a bounded 
measurable function on Q. Consider the functional 



on sets containing U and compactly contained in fi. Suppose E minimizes this func- 



1. If dU is C 1 , then BE is a C 1 submanifold ofQ. 

2. If dU is C 1 ' 01 , < a < 1/2, then dE is a C 1 '" submanifold of Q. The C x > a 
estimates depend only on the distance to dfl, esssup|/|, C 1,a bound for dU , 
and C l bounds on the metric in Q. 

3. If dU is C 2 and f = 0, then dE is C 1 ' 1 , and C°° where it doesn't touch U . 

Our initial value formulation falls into this category of problem. So our iV f 's and 
iV t + 's are C 1,a . Furthermore 



in local C 1 ^ convergence < f3 < a. 

The locations where N t ^ correspond to the jumps discussed earlier. To 
examine these areas we need to introduce minimizing hulls. 




(4.10) 



tional. 



lim N s = N t 



lim N s = N+ 



(4.11) 
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Definition 4.2.5. Let f2 be an open set. We call E a minimizing hull if 

\d* E r\ K\ < \d* F n K\ (4.12) 

for any F containing E with F \ E in K a compact set in Q. We say E strictly 
minimizes if equality implies that E and F agree in fl up to measure zero. 

The intersection of (strictly) minimizing hulls is a (strictly) minimizing hull. So, 
given a set E we can intersect all of the strictly minimizing hulls which contain 
E. This gives, up to measure zero, a unique set E' that we will call the strictly 
minimizing hull of E. Since E' is strictly minimizing, E" = E'. 

Solutions to the initial value problem given by equation ( 14.91) . have level sets that 
are minimizing hulls as follows. 

Lemma 4.2.6. Suppose that u is a solution to (|4.9p . Assume that M has no compact 
components. Then: 

• For t > 0, E t is a minimizing hull in M . 

• For t > 0, Ef is a strictly minimizing hull in M . 

• For t > 0, E' t = Ef if E^ is precompact. 

• For t > 0, \dE t \ = \dE^~\, provided that Ef is precompact. 

• Exactly when E is a minimizing hull \dE \ = IdE^l 

This is "Minimizing Hull Property (1.4)" in [7]. These minimizing properties 
characterize how the weak flow differs from the classical flow. The classical flow 
runs into trouble when the mean curvature changes sign. If the mean curvature is 
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negative on a part of N t , we could decrease the area of N t be flowing that patch out. 
So the weak flow avoids these areas by making sure that is always the outermost 
surface with its area. This sometimes necessitates jumping. 

The existence and uniqueness of solutions to equation (I4.4p with initial data are 
beyond our scope. However, for completeness, here is the existence and uniqueness 
theorem (3.1) from [7J. 

Theorem 4.2.7. Let M be a complete, connected Riemannian n-manifold without 
boundary. Suppose there exists a proper, locally Lipschitz, weak subsolution of (14. 9 p 
with a precompact initial condition. 

Then for any nonempty, precompact, smooth open set E in M, there exists a 
proper, locally Lipschitz solution u of (14. 9 p with initial condition E , which is unique 
in M\Eq. Furthermore, the gradient of u satisfies the estimate 



for each < r < a(x). 

The function a(x) depends on the Ricci curvature of M near Eq, but is always 
positive. More importantly, the requirement for a subsolution is satisfied by any 
asymptotically flat manifold. A function like ln(R) in the asymptotic end will suffice. 



Now that we have outlined the flow, we will describe some useful properties. 

Lemma 4.3.1. Let (E t ) t>0 solve (14. 9p with initial condition E . As long as E t 
remains precompact, we have the following: 




sup 

dE nB r (x) 




(4.13) 



r 



4.3 Useful Properties of Weak IMCF 
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• e * \dE t \ is constant for t > 0. 

• If Eq is its own minimizing hull then \dE t \ = e* |9i?o|- 

Proof. Since each E t minimizes the same functional, they must all have the same 

value for J u (E t ). Applying the co-area formula to the integral in J u gives 

J u {E t ) = \dE t \ - / -L / |V«| dAds (4.14) 
Jo |Vu| J dEs 

= \dE t \ - [ \dE s \ ds. (4.15) 
Jo 

Which has solutions of the form Ce l . By the minimizing hull properties, we could 
replace dE t with dEf for t > 0. Since dE^ is the limit of dE s as s \ 0, C = IdE^l. 
If Eq is its own minimizing hull then \dEo\ = \8Eq |, and C = \8Eq\. □ 

The next two lemmas tell us that when the classical solution exists, it agrees with 
the weak solution. 

Lemma 4.3.2. Let {N t ) c <t<d be a smooth family of surfaces of positive mean curva- 
ture that solve ( 14. ip classically. Let u = t on N t , u < c inside N c , and E t = {u < t} . 
Then for c < t < d, E t minimizes J u in Ed\E c . 

Lemma 4.3.3. Let Eq be a precompact open set in M such that OEq is smooth with 
H > and Eq = E' . Then any unique solution (E t )o < t<oo of (14. 9p with initial 
condition E coincides with the unique smooth classical solution for a short time, 
provided E t remains precompact for a short time. 

The authors of [7J point out that the stopping point for these theorems will be 
when either E t is no longer a minimizing hull, the mean curvature goes to zero, or 
the second fundamental form is unbounded. 
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The proof of existence and uniqueness of the weak flow are beyond the scope of 
this thesis. However, the method is as follows. First they assume the existence of a 
subsolution v. Then they solve the regularized problem 



IWI 2 + e 2 



E e u e := div | , ' = | - \l |W| 2 + e 2 = in tt L 



(4.16) 



u e = on 8E 

u e = L - 2 on dF L . 

Where Fl = {v < L}, e is small, and L is large, but bounded in size by a function 
of e. Then they take the limit as L — > oo, e — > 0. Some estimates on |Vw| and H 
guarantee that the solution passes to the limit. 

The authors also proved that in the asymptotic regime, the Hawking mass of the 
level sets of the flow converges to the ADM mass of the manifold 

Theorem 4.3.4. Assume that M is asymptotically flat and let (E t ) t > to be a family 
of precompact sets weakly solving (14.1 p in M. Then 

lim m H (N t ) < m ADM {M). (4.17) 

t— >oo 

Proof. They show that N t must approach coordinate spheres in the asymptotic 
regime. This is their Lemma 7.4. □ 

4.4 Geroch Monotonicity 

The Geroch Monotonicity formula says that the Hawking mass is nondecreasing 
under the inverse mean curvature flow. The original use of this was to propagate the 
mass of a surface out to infinity to compare to the ADM mass using IMCF. However, 

49 



it also provides bounds on integrals of H and area near the surface. We will first 
show how it arises in the smooth case, and then extend it over the jumps in the weak 
flow. 

In the smooth case, we simply recall that the Hawking mass is given by 



m H 



The first variation of area is given by 

^dA = Hr]dA. (4.19) 

Thus under smooth IMCF we have ^ gL4 = dA. The variation of H is given by 

dH 

— = A(-r]) -\A\ 2 r]-Rc{u,u)r]. (4.20) 

(JjL 

Thus if we look at the integral J H 2 dA under IMCF we get 

J H 2 dA = j H 2 - 2 |V ^f 1 + 2 \A\ 2 - 2 Rc(z/, v) dA. (4.21) 



dt 



The Gauss equation contracts to give 



K = K 12 + AxA 2 = | - Rc(i/, u) + \{H 2 - \A 2 \) (4.22) 



Here K and R are the scalar curvatures of iV and M respectively, K\ 2 is the sectional 
curvature of M in the plane tangent to N, and A« are the principal curvatures of N. 
Using this equation to cancel the Rc term gives us the following equation 

d r -t2 _ f J^nH\ 2 2 



dtlN H^] N 2K-2^- W -R (4.23) 
= JjK-R- 2^f- - \ (Ax + A 2 ) 2 - \ (X, - A 2 ) 2 . (4.24) 
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If R > 0, 



1 r H 2 = 4n X (N t ) - \ 



dt 



N 



H 



N 



N 



If iV* is connected, 



(4.25) 



d 
It 



H 2 < 8vr 1 



iV 



16tt 



N 



(4.26) 



In addition recall that \N t \ = \N \ e* in the smooth case. Thus we get that 

1 



/ d d 

VlQir—mn = — 
dt dt 



t/2l 1 _ 



167T 



H 



N 



(4.27) 



> 



V (l- — [ H 2 ]- 

2 V ^7lJ N J 167T 



1 



167T 



H 



N 



0. (4.28) 



Hence the Hawking mass is nondecreasing. 

To cover the gap, we simply note that at the jumps, the new surface E' t = Ef is 
a minimizing hull for E t . That means that where their boundaries differ, dE' t must 
have zero mean curvature (else a variation could decrease its area keeping it outside 
of Et, contradicting its minimizing property.) Hence 

H 2 > I H 2 . (4.29) 

dE t JdE' t 

Thus since \dE' t \ = \dE t \ for t > 0, we see that the, even at jumps, the Hawking 
mass can't decrease. This doesn't cover the possibility of dense jumping, or similar 
pathological behavior, but extending the Geroch formula to those cases requires using 
the elliptic regularization. 

The statement of Geroch Monotonicity given in [7] for one boundary component 
is as follows 
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Theorem 4.4.1. Let M be a complete 3-manifold, Eq a precompact open set with 
C 1 boundary satisfying j dEo \A\ 2 < oo, and (E t ) t>Q a solution to (14.91) with initial 
condition E . If E is a minimizing hull then 



m H (N s ) > m H {N r )- 



+ 



(16vr) 3 / 2 



/ 16tt - 8>Kx(N t ) + [ (2\D\ogH\ 2 + (\ 1 -\ 2 ) 2 + R) 

Jr L JN, 



dt (4.30) 



for < r < s provided E s is precompact. 
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Chapter 5 

Negative Point Mass Singularity Results 



5.1 Negative Point Mass Singularities and IMCF 

Although we do not need this fact, it is interesting to note that near a negative point 
mass singularity, we can define the inverse mean curvature flow. Even though there 
is no initial surface to start from we can take a limit of solutions to IMCF for starting 
surfaces that converge to p. This actually defines a unique solution u to the weak 
inverse mean curvature flow. This was shown in recent work by Jeffrey Streets [TO] . 

Theorem 5.1.1 (|10|). Let M 3 be an asymptotically flat manifold with finitely many 
singularities at {pi}. Then there is a unique solution to (14.411 on M \ {pi}. 

In this case since the area of the level sets is exponential in time, and surfaces 
near the singularities have vanishing area. Thus this flow only reaches the singularity 
at time — oo. 

Streets also showed that these surfaces were the best possible surfaces, in terms 
of the Hawking mass, as in the following theorem. 
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Theorem 5.1.2 ([10]). Let St be the family of hyper surf aces defining the solution to 
IMCF above. Let P t be any other family of hyper surf aces approaching the singularity. 
Then, 

lim m^(-Pi) < lim rriH(St). (5.1) 

t— >— oo t— oo 

We can also extend Geroch Monotonicity down to t = in the case where our 
initial surface has negative Hawking mass. 

Lemma 5.1.3. Let E be a surface in an asymptotically flat 3 manifold. Let E' be 
the boundary of the minimizing hull o/E. Let E or E' have negative Hawking mass. 
Then 



m#(E) < m H (T!). 



(5.2) 



Proof. If E' has nonnegative Hawking mass then toh(E') > > mn(S) and we are 
done. Thus we can assume that E' has negative Hawking mass. Since E' has negative 
Hawking mass, it must intersect E on a set of positive measure. Otherwise, E' would 



be a minimal surface, with Hawking mass y > 0. We define the following sets: 



E'\E Q 



s\s c 



(5.3) 



Recalling that |E + | < |E_| by the minimization property, and that H = on E 



+ • 



we observe the following: 



> 



E | + |E_ 

(16tt) 3 / 2 



igtt - / if 1 I u- } ///„(>:; 

So 



(5.4) 

(5.5) 

(5.6) 
□ 
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With this lemma and Geroch Monotonicity we can prove the following lemma. 

Lemma 5.1.4. Let (M,g) be an asymptotically flat manifold with ADM mass m, 
nonnegative scalar curvature and a single regular negative point mass singularity 
p. Let {Sj} be a smooth family of surfaces converging to p, which eventually have 
negative Hawking mass. Then for sufficiently large i, 



Proof. Since for large enough i, Ej has non-positive Hawking mass, we can apply 
Lemma r5.1.3l to show that E^ must have larger Hawking mass. From this surface, we 
start Inverse Mean Curvature Flow. Theorem 14.4.11 tells us that the Hawking mass 
of the surfaces N t defined by IMCF starting with E^ only increase. Theorem 14.3.41 
tells us that the increasing limit of the Hawking masses these surfaces is less than 
the ADM mass. Thus the Hawking mass of the starting surface was also less than 
the ADM mass. □ 

Now we relate the limit of the Hawking masses to the regular mass. 

Lemma 5.1.5. Let (M,g) be an asymptotically flat manifold with nonnegative scalar 
curvature and a single regular negative point mass singularity p. Then there is a 
smooth family of surfaces {Sj} converging to p such that 



m H (£i) < m. 



(5.7) 




(5.8) 



Proof. The Hawking mass of a surface Ej is given by 




(5.9) 
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Since the areas of the surfaces are converging to zero we have 



lim mH (S t ) = - lim -^EL f H 2 dA. (5.10) 

By the Holder inequality this is bounded as follows 

f H 2 dA<- ( l — ( [ H^dA]' 2 . (5.11) 



(16vr)3/2 y Si " (16vr)3/2 J 

Switching to the resolution space, we use the formula 

H = Jp~ 2 H + AJp- 3 V(lp). (5.12) 



Putting this into the previous equation we get 



/ H^ 3 dA= [_ (Tp- 2 H + 4^- 3 I7(^)) 4/ VZ4 
is, is, 



(5.13) 



= / (ipH + AV(Tp)) 4/3 dA. (5.14) 

Since Tp is zero on S and H is bounded, the first term goes to zero. The second term 
converges since the family of surfaces are converging smoothly. 

lim/ (TpH + 477(^)) 4/3 dA = 4 4 / 3 / ' Vi^Y^lA. (5.15) 



Combining all of these equations we have 

\ 3 / 2 



limm H (£i)<-~ (- [ V{Tp)^dA] 



m R (p). (5.16) 



To see when this estimate is sharp, we look at inequality (15. lip since that is the only 
inequality is our estimate. In the limit, this inequality is an equality exactly when 
the ratio of the maximum and minimum values of H approaches 1. We choose a 
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resolution such that f((p) = 1 on the boundary. We also choose a family of surfaces 
Sj given by level sets of Tp. Then if we look at the ratio 

lim ?™ = lim + ^ , (5.17) 

¥>->o H max </>-+o max + Au(ip) 

and remember that H is bounded, we see that the V{Tp) terms dominate, and as 
Tp — > 0, this ratio approaches 1. Thus with this resolution and this family of surfaces, 
inequality (15. 16[) will turn to an equality. □ 

With these results we can prove the following theorem 

Theorem 5.1.6. Let (M,g) be an asymptotically flat manifold with nonnegative 
scalar curvature and a single regular negative point mass singularity p. Then the 
ADM mass of M is at least the mass of p. 

Proof. First consider the case when p can be enclosed by a surface, S, with nonneg- 
ative Hawking mass. The minimizing hull of a surface with nonnegative Hawking 
mass has nonnegative Hawking mass. Thus we can run IMCF from £', and the 
AMD mass of M is at least mn(S') > 0. However, the regular mass of p is always 
nonpositive so in this case we are done. 

Now assume that p cannot be enclosed by a surface with nonnegative Hawking 
mass. By Lemma 15.1.41 we know that the ADM mass is greater than the Hawking 
masses of any sequence of surface converging to p which have negative Hawking 
mass. By Lemma 15.1.51 we know that there is a family of surfaces converging to p 
which have the mass of p as the limit of their Hawking mass, hence the ADM mass 
is greater then their Hawking masses which limit to the regular mass. □ 
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This can be extended to general negative point mass singularities. However, first 
we need to consider the effect of multiplication by a harmonic conformal factor on 
the ADM mass of a manifold. 

Lemma 5.1.7. Let (M 3 ,g) be an asymptotically flat manifold. Let be a harmonic 
function with respect to g with asymptotic expansion 

<p=l + ^- + 0\ | . (5.18) 

Then, if the ADM mass of (M 3 , g) is m, the ADM mass of (M 3 , (p 4 g) is m + 2C . 

Proof. This is a direct calculation. We write g v = gip A , and calculate, only keeping 
the terms of lowest order in 1/ \x\ since we are taking limits as \x\ — > oo. 

m v = , l } m 77T I (9ij,i ~ 9u,j) n J dA (5.19) 

\x\-Hx l07T J s s J J 



\x 



¥ A f <P* f 

lim — - / (//,,., - </„.;) n> (I A ■ lim — / - 5 u <pj) n J dA (5.20) 

5 1 -+oo 107T J s s ' |a;|->-oo 47T J $8 

lim ip 4 m + lim <p 3 lim — / (ipj — 3ipj) n? dA (5-21) 

»oo \x\— >oo \x\— ¥00 All J gs 

m+ lim — / (pj-nPdA (5.22) 

m+ lim — / (V(p,v) dA (5.23) 

m + 2C. (5.24) 



□ 



Using this we can now extend Theorem 15.1.61 to a general singularity. 
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Theorem 5.1.8. Let (M,g) be an asymptotically flat manifold with nonnegative 
scalar curvature and a single negative point mass singularity p. Then m, the ADM 
mass of M, is at least the mass of p. 

Proof. If the capacity of p is nonzero, then the statement is trivial. Thus we assume 
the capacity of p is zero. Using the terminology of Definition 12.3.41 Theorem 15.1.61 
tells us that the ADM mass of (M, hfg) is at least the mass of the regular singularity 
at Ej = Pi. Each hi is defined by the equations 

Ah,, = (5.25) 

lim hi = 1 (5.26) 

x— >oo 

hi = on (5.27) 

Thus it has asymptotic expansion 

hi = l-Q + (-^j . (5.28) 

Where Cj is the capacity of Ej. Thus, the ADM mass, rrii, of (M, hfg) is given by 
m — 2Cj. Now we know that m 8 > m^pi). Taking lim of both sides gives us 

lim rrii > lim m^(pi) (5.29) 

i— >oo i— >oo 

Since Cj is going to zero, the left hand side is simply m, and so has no dependence 
on which {Ej} we chose in our mass calculation. Thus we get 

m > sup lim m R (pj). (5.30) 

{Si}*-* 00 

as desired. □ 
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5.2 Capacity Theorem 

Perhaps the most important new results of this thesis are Theorems 15.2.41 and 15.2.51 
which relate the capacity of a point to the Hawking masses of surfaces near that 
point. 

The capacity of a surface provides a measure of its size as seen from infinity. We 
extend the definition of the capacity of surface to the capacity of a negative point 
mass singularity. We then show that if a NPMS has non-zero capacity the Hawking 
mass of any family of surfaces converging to it must go to negative infinity. First the 
definition of capacity: 

Definition 5.2.1. Let S be surface in an asymptotically flat manifold M. Define 
the capacity of £ by 



L J m ) 

It is worth noting that if £ and £' are two surfaces in M so that S divides M 
into two components, one containing infinity and the other containing £', then 



since the infimum is over a larger set of functions. 

We will next define the capacity of a singular point. The natural definition is the 
one we want. 

Definition 5.2.2. Let p be singular point in an asymptotically flat manifold M. 
Chose a sequence of surfaces S, of decreasing diameter enclosing p. Then define the 
capacity of p by the limit of the capacities of Sj. 




(5.31) 




(5.32) 
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Before using this definition we have to show that it is unique. 

Lemma 5.2.3. Let Ej and Ej be two sequences of surfaces approaching the point p. 
//limC(Ej) = K, limC(Ej) = K. Hence C(p) is well defined. 

Proof. Since the Ej are going to p, for any given Ej , we can choose Iq such that for 
all i > Zq, Ej is contained within E io . Thus if ip is a capacity test function for E io , 
i.e. y2(Ej ) = 1 and ip — > at infinity, then ip is also a capacity test function for Ej. 
Since C(Ej) is taking the infimum over a larger set of test functions than C(E io ), 
C(Ej) < C(E io ). Thus if we create a new sequence of surfaces Ej, alternately choosing 
from Ej and Ej, such that each surface contains the next we get a nonincreasing 
sequence of capacities. Thus if either original sequence of surfaces has a limit of 
capacity, then this new sequence must as well, and it must be the same. Hence, 
linij^oo C(Ej) = linij^oo C(Ej). □ 

Now we look at the relationship between capacity and the Hawking mass of a 
surface. We will use techniques similar to those used in PJ. 

Theorem 5.2.4. Let M be an asymptotically flat 3 manifold with nonnegative scalar 
curvature, and negative point mass singularity p. Let Ej be a family of surfaces 
converging to p. Assume each Ej is a minimizing hull. Assume the areas of Ej 
are going to zero. Then if the Hawking mass of the surfaces is bounded below, the 
capacities of surfaces of foliation near p must go to zero. 

Proof. To use Geroch monotonicity, we need to know that our IMCF surfaces stay 
connected. In the weak formulation of IMCF, the level sets E t always bound a 
region in M. Thus if E t is not connected, one of its components EJ 1 must not bound 
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a region. That is, E£ is not nomotopic to a point in M. Since M is smooth, it 
must have finite topology on any bounded set. Thus we know that near p, there is a 
minimum size for a surface that does not bound a region. Call this size v4 min . Thus 
if we have any surface that does not bound a region, it must have area greater then 
A m i Q . According to Lemma T4.3. II the area of our surfaces grow exponentially. Thus if 
we restrict ourselves to starting IMCF with a surface with area A min /e, and only run 
the flow for time 1, we will stay connected. At first glance it seems we may need to 
worry about the jumps in weak IMCF, however Geroch monotonicity doesn't depend 
on smoothness of the flow, and neither does the area growth formula. Lemma 14.3.11 
holds from t = 0. Thus even with jumps, the area of our surfaces will remain below 



Here, the integral is only over the portion of M outside of E. Call this integral, £(p), 
the energy of p. Thus for any (p with p(oo) = and <£>(E) = 1 we have £{p) > C(E). 
So we will find an estimate that relates the Hawking mass and the energy of a test 
function p. 

Choose a starting surface E with sufficiently small starting area. Let / be the 
level set function of the associated weak IMCF starting with the surface E. Call the 
resulting level sets E 4 . Now if we use a test function of the form <p = u(f ), then the 
energy of p is given by 



.4 



mm • 



Now recall that capacity of a surface is defined by 




(5.33) 




(5.34) 
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Since / is given by IMCF, we know that ||V/|| = H where H is the mean curvature 
of the level sets. Next we will use the co-area formula. This states that if we have 
a function zona domain fl, and a function h : R — > R so that the range of h(z) is 
[a, b], then 

I hdV = [ h(z(t)) [ \Vz(t)\ dA t dt. (5.35) 

Jn Ja J St 

Here S t are the level sets of h(z(t)). If we use the co-area formula with the foliation 
S t , our integral becomes 

E(<p)= [ (u'(t)) 2 I \H\ dA t dt. (5.36) 

Here the co-area gradient term cancels one of the ||V/|| terms. Now we will bound 
the interior integral of curvature. We know that IMCF causes the Hawking mass to 
be nondecreasing in t. We first rewrite the definition of the Hawking mass mn(SJ) = 
m(t) as: 

J H 2 dA t = 16tt ^1 - m(t) . (5.37) 

Here A(t) is the area of T, t - Since the Hawking mass is nondecreasing under IMCF 
we have: 

j H 2 dA t < 16tt ^1 - m(0) y^^) • (5-38) 



Thus we can use Cauchy-Schwartz to get: 



J \H\dA t <^W)^ 167T ^l-m(0)^^j. 



(5.39) 
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We can rewrite this as: 



J \H\ dA t < \JaA{t) + Py/A(j). (5.40) 



Furthermore, since A(t) grows exponentially in t, we can write this as: 

J \H\ dAt < \fae t + /3e*/ 2 = v(t). (5.41) 

Where A has been absorbed into a and j3. Thus our energy formula has become 

roc 

£(<p) < / (u'(t)) 2 v{t) dt. (5.42) 
Jo 

with 

v(t) = yfae* + fie 1 ' 2 (5.43) 

where a = IQtxAq, f3 = (167r) 3 / 2 A^j 2 \m \, and A is A(E Q ). This means we can pick 

our test function u(t) to be as simple as: 

. 1 -f < t < 1 

u(t) = t ~ ~ (5.44) 



t > 1 



Then our integral becomes: 

E(ip) < [ v(t)dt 
Jo 



\fae l + /3e*/ 2 dt 

= [ e^Vae'/ 2 + (3dt 
Jo 

re 1 /* 

= 4 J a/ ax 2 + (3 dx 

re 1 /* 

< 4 J ^x + 



= 2V^(e 1/2 -l)+4 v ^(e 1 / 4 -l) 
< 2^ + 2^. 
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Since mn(S) < y j^, m is bounded above. By assumption m is bounded below, 

so a and (3 are bounded by multiples of A and y/A~o respectively. Thus £(<p) goes to 
zero if Aq — > and mo is bounded. Hence C(p) must be zero since it is the infimum 
over a positive set with elements approaching zero. □ 

Theorem 5.2.5 (Capacity Theorem). Let M be an asymptotically flat 3 manifold 
with nonnegative scalar curvature, and negative point mass singularity p, such that 
there exists a family of surfaces, E i; converging to p with area going to zero. Then if 
the capacity of p is nonzero, the Hawking masses of the surfaces Ej must go to — oo. 

Proof. Any such family of surfaces will generate a family of minimizing hulls that 
will also converge to p. By Lemma [5.2.4[ the masses of {S'j} must go to — oo. Thus 
the masses of {S^} must go to — oo. Thus for sufficiently large i, the masses of the 
minimizing hulls are all negative. From then on Lemma [5.1.31 applies, and the masses 
of Ej must be less then the masses of E^. Hence they also converge to — oo. □ 
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Chapter 6 

Symmetric Singularities 



In this chapter we will look at some more examples of negative point mass singu- 
larities. The symmetry ansatz provides more structure than in the definition of a 
NPMS. 

6.1 Spherical Solutions 

A spherically symmetric manifold, M, has a metric given by 

ds 2 = dr 2 + ^P-dS 2 . (6.1) 

We can directly calculate that the scalar curvature of M is given by 

„ 16nA + A' 2 - AAA" . , 

R= ^ • (6.2) 

For this manifold to be asymptotically flat it is necessary for the Hawking masses of 
the coordinate spheres to approach a constant. The Hawking mass of a coordinate 
sphere is given by 
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Due to the spherical symmetry of the manifold, we know that if IMCF is started 
with coordinate spheres, it must continue with coordinate spheres. Thus we know 
that this quantity must be non- decreasing. For this manifold to be asymptotically 
flat, this quantity must have a limit at oo. This limit is the ADM mass. 

The only possible location for a singularity in such a manifold is at the origin. 
We can find a straight forward function to resolve the singularity. We need a smooth 
function such that 

lim = 4vr, (6.4) 

r— ¥00 {p 4 

or any other constant. 

Multiplying by ip~ A to find the model space changes our metric to the form 

ds 2 = dp 2 + ^W. (6.5) 
Atx 

Where A goes to 4ir as p approaches zero. Let £ be the surface p = 0. The behavior 
of ip away from the singularity is very flexible as long it is bounded, nonzero, and 
goes to one at infinity. In order for this to be a regular singularity, we need that p 

be well defined. Thus we must require 

f r f dr 

P (r) = J q <p- 2 dr = ^ < 00 (6.6) 

for finite r. Any spherically symmetric singularity with this condition on A{r) must 
be regular. The regular mass is given by 



Since we have only defined the asymptotic behavior of ip near r = 0, we will compute 
everything as limits as r goes to zero. First we need to find v(ip). Since p = J Q r <p~ 2 dr, 
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a chain rule calculation of ^ yields 



2 dip 



u(<p) = <p 2 -£. (6.8) 



Putting this in and expanding the other terms in the definition of the regular mass 
gives us 

m R = -limi(i/ (---) -V,4].^ (CO) 



>o 4 



Using equation (16.41) and l'Hopital's rule we find that 

<9y? 1 <9A 
dr 167T(^ 3 <9r 

Continuing from above 



(6.11) 



r->-0 4 7T ' 25b7l z 



l im L^1^!a 3 / 2 A /2 (6.12) 

r^0l0247T 7 / 2 A 2 ^ ^ 



1 A' 2 

hm 



r ^ 647T 3 / 2 A 1 / 2 ' 



This agrees with the limit of the Hawking masses of coordinate spheres as r — > and 
hence A — )■ as well. 

For completeness, we can also examine the capacity of the central point in these 
solutions. The capacity of a coordinate sphere in an asymptotically flat spherically 
symmetric manifold is given by the 1/p term in of the harmonic function that is 1 
on the sphere and at infinity. Since harmonic functions have constant flux with 
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respect to the outward normal, we could also describe this function as the constant 
flux function which goes to at infinity and 1 on the coordinate sphere. Then the 
capacity of the sphere is given by the flux constant of this function, divided by — An. 
Reversing this definition we define the following function: 

Definition 6.1.1. Let (M,g) be an asymptotically flat spherically symmetric man- 
ifold. Let / be the radial function that has constant outward flux — An through 
coordinate spheres and goes to zero at infinity. Call / the radial capacity function 
for (M, g). 

This definition allows the following lemma. 

Lemma 6.1.2. Let f(r) be the radial capacity function for the manifold (M,g). 
Then the capacity of the coordinate sphere at r = r is given by /(r ) _1 . 

Proof. Consider the function f(r)/f(r ). This function is 1 on the coordinate sphere 
r = r , goes to zero at infinity, and is harmonic. In the asymptotic regime its l/p 
coefficient is l//(r ). Thus this is the capacity of the sphere r = r . □ 

Thus the capacity of the central point is given by lim r ^ /(r) -1 . To calculate this 
value we first note that 



df An 



(6.13) 



dr A(r) 
since / has constant flux — An. Thus 

/OO /-00 7 

f'(r)dr = An^ — . (6.14) 

In order for the capacity of the central point to be nonzero, this must be finite. 
Asymptotic flatness tells us that the part of the integral in the asymptotic regime is 
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finite. Thus the only concern is where r — > and hence A(r) — > 0. Thus if 



is finite, our central point has positive capacity. For example, if we assume that as 
r — 7- 0, A(r) is asymptotically a multiple of a power of r, as kr p . Then the capacity 
is positive exactly when p < 1. In this case, we see by equation ( 16.121) that the mass 
of the singularity is infinite, confirming Theorem 15.2.51 However ifl<r<4/3, we 
see that the mass of our singularity is still infinite, but the capacity is now finite. 
This removes the possibility of strengthening Theorem 15.2.51 into an if and only if 
without additional hypotheses. 

6.2 Overview of Weyl Solution 

We will be looking at axisymmetric static vacuum spacetimes. These examples have 
two Killing fields, one spacelike and one timelike. These reflect the rotational and 
time translation symmetry of our spacetime. Furthermore the timelike Killing field 
is hypersurface orthogonal, and the two Killing fields commute. Since our two vector 
fields commute, we can call the timelike field dt and the spacelike field dg. Since 
they are Killing fields, we know the metric is only a function of the remaining two 
coordinates Xi,x 2 . We can also assume that do is orthogonal to d t . Thus our metric 
is of the form 



Here A,B,g^ are functions of X\ and xi- We set x\ = AB, and chose xi to be 
orthogonal to the other coordinates. Renaming x\, p and X2, z, our metric takes the 




(6.15) 



g = —A 2 dt 2 + B 2 d9 2 + gndxl + g\idx\dx2 + gndx\. 



(6.16) 
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form 

g = -A 2 dt 2 + p 2 A- 2 d6 2 + U 2 dp 2 + V 2 dz 2 . (6.17) 

See Theorem 7.1.1 in [12] for an explanation of the lack of cross terms and further 
details on this derivation. 

We may also rescale z to set V = U . Our metric is now encoded in the functions 
A(p, z) and U(p, z). If we define A and p by 

A = In A p = \n(AU) (6.18) 

our metric looks like 

g = -e 2X dt 2 + e~ 2X [p 2 d6 2 + e 2 ^ (dp 2 + dz 2 )] . (6.19) 

Now, if we compute the curvature of this metric, and apply the vacuum condition 
we get the following equations for A and p: 

= A pp + -\ p + \ zz (6.20) 

p p = p (A* - Aj (6.21) 
p z = 2pX p X z . (6.22) 

The first one is the same as the statement that AA = when viewed as a function 
of flat M 3 with coordinates (r,z,9). We can use this to generate spacetimes. We 
can think of A as a classical potential function. However, we should not think that 
g gives the metric with this gravitational potential. For example, the Schwarzschild 
solution is generated by a A that is not spherically symmetric 

As always we are interested in asymptotically flat spacetimes. In this case we 
want our metric to be the flat metric in cylindrical coordinates at infinity. Thus we 
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need the following asymptotics on A and p: 

lim A = (6.23) 



r— >oo 



lim n = (6.24) 

Where r 2 = p 2 + z 2 . Since A is flat-harmonic, we know that at infinity it looks like 
C |r| _1 + 0{r~ 2 ). This decay and equation (I6.20p tells us that 

lim |A P |,|A Z | < %. (6.25) 

r— >oo r 

Putting those into equations (I6.2ip and (I6.22p gives us 

C 

lim \p p \ , \p z \ < — . (6.26) 

These conditions are enough for asymptotic flatness. Thus all that is required of our 
metric metric for it to be asymptotically flat is that A and \i approach zero at oo. To 
calculate the ADM mass of g, we need coordinates that are asymptotically Cartesian 
rather then asymptotically cylindrical. Using the change of coordinates 

x = p cos 6 y = psmO (6.27) 

we get that our metric is 



2A , 2 2 a A (e 2M - 1) x 2 \ , 2 e~ 2A (e 2 ^ - 1) xy , , 

g = - e 2X dt 2 + e~ 2X 1 + ^-5 dx 2 + K - ^-dx dy 

\ x z + y z J x z + y z 

+ e -2A ( ! + (e 2 ^ -l) y 2 A 2 + e -2A e 2 M ^ 2 (6 2g) 
\ x 2 + y 2 J 

If we plug this into the formula for the ADM mass we get the following: 

r e -2A 

m = lim / [1 - 2zA 2 - 2pA„ + e 2 ^ (2pu p + 2zp z - 2p\ p - 2z\ z - 1)1 dA s 

(6.29) 
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As r grows the terms X z and X p are at most order r 2 . The derivatives of p, are at 
most order r -3 . The function e~ 2A is 1, as is the function e 2fl . Therefore this integral 



is 



m = — lim I -(z\ z + p\ p ) d A s = lim [ (u,VX) 5 dA 5 . (6.30) 

All r— >oo J g T All r^foo J g 

Since A is harmonic in the flat metric, we can compute this on any surface homotopic 
to a large sphere at infinity. 

Lemma 6.2.1. Let (M,g) be the t = slice of an asymptotically flat axisymmetric 
vacuum static manifold with metric 

g = e ~ 2X [p 2 d6 2 + e 2 ^ {dp 2 + dz 2 )] . (6.31) 

Then the ADM mass of (M, g) is given by 

~ lim / (is,V\) s dA s , (6.32) 

47T r-5>oo J 

where E is any surface enclosing the singularities of X. 

As in the spherically symmetric useful harmonic function can tell us about 

the capacity of the singular points. In this case our metric comes with a harmonic 
function. The function e A is harmonic in our metric. It goes to 1 at infinity and 
infinity at any positive singularities of our potential. Thus the capacity of level sets 
of e A must go to zero as A — > oo. This tells us that these singularities have zero 
capacity. 

6.3 Zippoy— Voorhees Metrics 
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The particular family of metrics we will consider are the Zipoy-Voorhees, or 7 
metrics. These are given by the potential arising from a uniform density rod of length 

2a and mass m at the origin. This gives a potential of 

m R + + R„- 2a m 2 AR + R_ 

X= 2-a ln R + + R. + 2a » = (R + + _R_) 2 — 4a 2 ' (6 - 33) 



Where R± = p 2 + (z ± a) 2 . If m = 2a, then this metric is Schwarzschild. It 
only represents the area outside the horizon. The interval [—a, a] on the z-axis is 
the event horizon. If m 7^ 2a then this metric has a naked singularity at p = 0, 
\z\ < a. See [Ej for more information on the cases with positive m. When m^O, 2a, 
the resulting spacetime has ADM mass m and the [—a, a] interval on the z-axis 
becomes a candidate for a negative mass singularity. To investigate the area near 
the singularity, we will look at p constant cylinders from z = —a to a. The area of 
these cylinders is given by 

2np I e^ 2X dz. (6.34) 

J —a 

Since both A and /1 are given by logs, the integral simplifies to 

W-A (fl> + «-)■-*»■ J U + fi- + a. J iz - (635) 

When z is between —a and a, and p is small we have the following approximations: 

iR + R - ( « 2 -- )2l +0 (l) (6.36) 



(R + + R_) 2 -4a 2 a 2 p 2 

-p 2 + 0( P <). (6.37) 



R+ + R--2a 1 2 , 



i?+ + R- + 2a 4(a 2 - z 2 ) 
Hence our integral becomes 
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Now if m 7^ a, the p term has positive exponent. Thus the areas of these surfaces 
go to zero. The fact that the function e A is harmonic and goes to infinity near the 
[—a, a] on the z axis tells us that these surfaces have zero capacity. Thus, as long as 
m ytz a, these fulfill the definition of a negative point mass singularity. 

Continuing in this fashion we can estimate the mass of this singularity. We will 
just estimate the mass of the singularities using the level sets of the function A. 
Referring back to Definition 12.3.41 our function hi given by the level set A = L is 
given by 

hi = — — . (6.39) 

L- 1 L- 1 v ; 

Here L = ±i, with the sign chosen to be the opposite sign to m. We can calculate 

z/(^) 4 / 3 as 



uW = e l*"fr (a; + Xlf 3 . (6.40) 

(L — 1) 



Now we will approximate the surface A = L by a level set p = pi. Noting that A p 
is much larger then X z tells us that this assumption is valid. With that assumption, 
our mass integral becomes 

£ = 7T^W e " fA f e ' ¥ K + ^) 2/3 ^. (6.41) 

(6 — 1 J J —a 

Now we note the first order behavior of A and p near p — 0, \z\ < a. 

,2 



m m 
A~ mp a ~ — — In p 

a ' a 1 



ml mz 
A„ ~ A, 



(6.42) 



a p a(z 2 — a 2 



Using these we approximate the above integral 
2vrp 

(jj-m/a _ 1) 



2np 2 m f a im* fm 2 2 m 2 z 2 A 2/3 , 
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As p goes to zero, this has p dependence as 



o 

— m 1 221 1 

£ ~ c ■ — m- ( 6 - 44 ) 

(p-m/a _ X ) 4 /3 



For some constant C. If m > 0, then the bottom term contributes a p»~ to the 
growth giving an overall power of |^ + ^ — 1. Otherwise it contributes negligibly. 
Thus we have the following p dependence 

C + p3^+ a m>0 , 

(6-45) 

C_p3^" 3 a m < 0. 



The exponent on the p is negative when ^ G ^— 1, — |J and negative when 
^ falls outside the closure of that range. Outside that range we have produced 
an example of a set of surfaces which give zero to the mass under Definition 12.3.41 
Hence, since that mass is a sup over all such surfaces, we know it must be zero. For 
^ inside that range, our surfaces give a mass of — oo. This is inconclusive, since it is 
entirely possible that there exists a better behaved family of surfaces. 

Checking our two known cases, — = ±1 we see that the positive Schwarzschild 
metric doesn't have a singularity, and the mass of the negative Schwarzschild is 
nonzero and finite. It shouldn't be surprising that for positive m outside that range 
we do not get a negative mass as a Negative Point Mass Singularity since these 
singularities are the only points without zero scalar curvature in a static manifold 
with positive ADM mass. It seems sensible that they shouldn't be assigned a negative 
mass. 
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Chapter 7 

Open Questions 



There are a number of unanswered questions having to do with negative point mass 
singularities. The most prominent is extending Theorem 15.1.81 to include multiple 
singularities. Bray, in [6], has a solution that depends on an unproven geometric 
conjecture. A further result would be what I have been calling the "Mixed Penrose 
Inequality." This would be a result that combines the two cases, singularities and 
horizons, and provides a lower bound on the ADM mass of a manifold containing 
horizons and singularities. The desired inequality is presented in Appendix [Bl 

It is clear that a removable singularity should have mass zero. Precisely what con- 
ditions on an negative point mass singularity guarantee that it is removable requires 
further investigation. 

On a more concrete note the Zipoy-Voorhees metrics with ^ e ^— 1, — |j 
deserve further investigation. In particular the metrics with ^ in that range and 
positive seem to offer the possibility of negative mass singularities in a vacuum static 
manifold with positive ADM mass. 
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On a broader scale, the current presentation of negative point mass singularities 
is only in the Riemannian case. Extending the definitions to the full Lorentzian 
context would require some equation to control the behavior of the manifold in 
the neighborhood of the singularity in the timelike direction. Furthermore, a new 
definition of mass would have to be added, since the current definition depends on the 
foliation of surfaces converging to the singularity. Furthermore, the conformal factor 
definition of the regular case doesn't translate in an obvious way to the Lorentzian 
case even with the Schwarzschild metric. Studying negative point mass singularities 
in a spacetime is an important direction to pursue. 
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Appendix A 



Miscellaneous Calculations 



A.l Calculation of total magnification of a NPMS lens 



The magnification of a negative point mass singularity lens is given by 

1 II l|4 2 
H=- ^T= jj 4 i =1 + ^ 2' ^ 



The image locations x± associated to a given y are 

x± = X - (y ± vV + 4m) . (A.2) 



First we calculate ||a;±|| 4 as 



|:r ± || 4 = i; (</ ± v / l/ 2 + 4m) 4 (A.3) 

= (y 2 ± 2f/vV + 4m + y 2 + 4m) (A.4) 

= — f4y 4 + 16y 2 m + 16m 2 + 4y 2 (y 2 + 4m) ± (2y 2 + 4m)4?/vV + 4m) 
16 V / 

(A.5) 

= (% 4 + 32y 2 m + 16m 2 ± (% 3 + 16ym) vV + 4m) (A.6) 

= ^y 4 + 2y 2 m + m 2 ± Qy 3 + ym^j \/ y 2 + 4m. (A.7) 
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Now we define the following terms 

1 



A=^y 4 + 2y 2 m (A.8) 
B = m 2 (A.9) 

C = y Qy 2 + mj vV + 4m. (A.10) 



Thus ||a;±|| 4 = A + B ± C. The negative image has negative magnification, so we 
have to look at [J,(x+) — fx(x-). This is given by 



*(*)= „ "l + " 4 2 - „ 2 (A.11) 

\\x + \\ — m A IF-!! — m z 



A + B + C _ A + B-C 
A+C A-C 



-2BC 
A 2 -C 2 ' 



(A.12) 



- l + ^C~i l + ^ (A.13) 



(A. 14) 



We can simplify this as 

2BC = ym 2 (y 2 + 2m) vV + 4m (A. 15) 

A 2 - C 2 = Qy 4 + 2y 2 m^j - y 2 ( Qy 2 + rn) vV + 4m) (A. 16) 

= -y 4 m 2 - 4y 2 m 3 = -y 2 m 2 (y 2 + 4m). (A. 17) 

Thus we get 

„ = 4±|^. (A.8) 
?/A/r + 4m 
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A. 2 Solutions to Cusp Equation 

In the following we will drop the * subscripts on the quantities Z, 7, and m. We are 
solving the equation 

= Z + e K (^ + ^jZ, (A.19) 

where 

Z = -4* (7 + -) =3, (A.20) 

and 



z = ±J — . (A.2i; 



Thus we get the following 



m • m 

~ = < ' * 7 + =2 

2^ z 



m3/2 - (e- - 7 ) 3/2 ^ = = - (e* - 7 ) 3/2 . 



(A.22) 



Z 3 ' z" z" 



Note the conjugate bar on the m on the last relation. Since m is negative m 3 ^ 2 
is imaginary, and so we had to multiply by —1 when we conjugated the previous 
equation. Thus our equation is 

Multiplying by 3^ and using the substitutions in ()A.22j) and then multiplying by 



e %Lp l 2 gives us 



= e-v* (e^ - 7) 3/2 + e K e^ (e^ - 7)^ . (A.24) 



If we define w = e %ip [e tLp — 7). Then we have 

= w 3 / 2 + e K w^ 2 . (A.25) 
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Which is solved when u> 3//2 purely real or imaginary, when e K is negative or positive 
respectively. This corresponds to w 3 being purely real and having the opposite sign 
as e K . The imaginary part of w 3 is 

Im(w; 3 ) = 7 sin cp [47 2 cos 2 tp — 67 cos cp + 4 — 7 2 ] . (A. 26) 

Setting aside the case when 7 = 0, we are left with the roots 



<Pi = 0, (f 2 = it, (p 3 A = acos , tp 5 = 2n - cp 3 , cp 6 = 2rr - ip 4 . 



47* 



(A.27) 



If 7 2 < 3/4 we only have at most 71 and 72. Looking at the real part of w 3 now 

w 3 ) = -4 7 3 cos 3 (p + 67 2 cos 2 ip + (37 3 - 37) cos <p - + 1. (A.28) 



Thus since there are no sin ip terms, the real part of w 3 only depends on cos (p. For 
each (pi, we only have a solution, and hence a cusp, for one value of e K , where e K has 
the same sign as w 3 . Thus we need to find where the sign of w 3 changes. All this is 
in Table [O 





COs(v9,) 


Re(V) 


+ 




<P1 


1 


-7 3 + 3 7 2 - 3 7 + 1 = (1 - 7 ) 3 


7 < 1 


7 > 1 




-1 


7 3 + 3 7 2 + 3 7 + 1 = (1 + 7 ) 3 


always 




^3,5 


3+ a/47 2 -3 
47 


1 ~ l7 2 + kl V47 2 - 3 


7 > 1 


V / 3?4<7< 1 


^4,6 


3+ a/47 2 +3 
47 


1 - |7 2 - \l V47 2 - 3 




always 



Table A.l: Values of the real part of w 3 for various 93 j. 
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Appendix B 



Work Towards the Mixed Penrose 
Inequality 



B.l Purpose 

This appendix is a record of my attempts towards solving what we have been calling 
the "Mixed Penrose Inequality." The desired conjecture is as follows: 

Conjecture B.l.l (Mixed Penrose Inequality). Let (M,g) be an asymptotically flat 
manifold with outer minimizing boundary £ and negative mass singularities p i . Also 
assume that (M, g) has nonnegative scalar curvature. Then 



The exponent on the negative point mass singularity term comes from consoli- 




(B.l) 
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dating the integrals for the masses as follows (assuming they are regular): 



(\ / " 3/2" 

X>N 2/ j =-lE[jg/ Ei ^) 4/ ^) 



2/3 N 



3/2 



(B.2) 



3/2 

11 



L"M i/3 da) (B.3) 



4 Vtt 



7T 

\ 3/2 

z%) 4/3 da . (B.4) 



This has so far been out of our reach. As a first case we have been working on 
a reduced conjecture which assumes the singularity and black hole have equal and 
opposite masses. Then all that is expected is that the ADM mass is positive. 

Conjecture B.1.2. Let (M,g) be an asymptotically flat 3-manifold, with minimal 
boundary £ and negative mass singularity p of opposite mass. Also assume that 
(M,g) has nonnegative scalar curvature. Assume that p can be resolved by a har- 
monic conformal factor ip so that (M, g) is asymptotically flat, £ is still minimal, 
and U, the resolution of p, is minimal as well. We can assume that ip = 1 at infinity, 
ip = only on U, ^(<£>) = 1 on H, and v{(p) = on E. 
Then the ADM mass of (M, g) is non-negative. 

There are three methods we have for proving Penrose style theorems. Inverse 
mean curvature flow, as in [7], which only takes into account a single horizon or 
singularity, so probably won't be helpful here. Bray's conformal flow of metrics as in 
[5] is the method we worked most with. Lastly, Bray's Minimal Surface Resolution 
technique used in [6], might also be useful, however to use it, we would need to adapt 
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it to accommodate the black holes. Furthermore, one of the steps in its proof is not 
yet complete. 



B.2 Effects of Harmonic Conformal Flow 

The basic idea here is to use a harmonic conformal flow on the model space (M,g), 
which is equivalent to a harmonic conformal flow on the actual space (M,g). 

If we look at a flow on ip, and set ip = ^<p, then we can look at the changes to 
the various quantities. If we set m to be the ADM mass, M to be the mass of the 
black hole, and N to the mass of the negative mass singularity. 



3/2 



M(E) 



16n 



1 - 



167T ./yr 



H 2 dA 



(B.5) 



(B.6) 



For a first estimate we will look at the rate of change at time zero for initial vp. First 



we look at II. A direct calculation shows 



a/3 



dA . (B.7) 



Under the assumption that we are working at t — 0, we can put in the known 
information about <p to get: 

~4 



^ ( n) = -?(i/ n K^) I/2 (i/ /M 



1/3 



-u(ip) + Hip 



dA 



N(U) 



J n 2v{^) + \Ej)dA 
IndA 



(B.8) 
(B.9) 
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If II is minimal as well, we can drop the H term. We can also drop this term if we 
assume that ip — on II. Likewise we can compute the change to M as 



1 




)""•(/, 



M(E) = V A dA) ( / 4>y 3 dA) (B.IO) 



= M(E) l^' (B11) 

The integral of H 2 term drops out since H is zero. We are also concerned with the 
change in the ADM mass of the manifold. We can measure that with the integral of 
the flux across a sphere at infinity, since this will measure the 1/r term of ip. (Here 
we assume that ip = at infinity.) Since the harmonic functions have constant flux, 
we know that the flux across S and II will also give us the 1/r term for if). Thus the 
change in m is given by 

rh= V(V>) -vdA + J V(^) -udAj (B.12) 

= ^(y ' vty)dA + J v{i))dA^ (B.13) 

Adding these together, the quantity we want to decrease is the mass surplus X, 

X = m- M - N (BAA) 

X = m- M - N (B.15) 
= i (/ u ^) dA + J v(i>)dA^) - (B.16) 

V kv'dA J V fn dA J 

If we can flow so that this decreases until iV (or M) is zero, then we can apply the 
Penrose (NPMS) inequality to get that X > at this time, hence X > at t = 0, 
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and we are done. At this point we tried different boundary conditions to accomplish 
this. 

B.2.1 Boundary Conditions 1 

For a first try, set ip by the boundary conditions i/j — at infinity, ip = — f y-^*r on 
£, with both integrals over S, and ip — and = — | on IT. This gives a value 
for X of 

X = — [ v{i))dA - — |n| - /3M - aiV. (B.18) 
B.2.2 Boundary Conditions 2 

As a second try, we can look at the boundary conditions, ip — at infinity, — a 
on S, //(■?/>) = /3 on IT, where we choose a to be — 2 j j.^ 3 , and pick /3 so that rh is zero. 
This means that the flux of ip through the two surfaces is zero. 

B.2.3 Unknown Boundary Conditions 

Trying the boundary conditions if) — on IT and at infinity while we leave the 
boundary conditions of ip on E up in the air for now, except that we want ip < 
and v(if)) > 0. If we set S to be a large sphere at infinity, then m will be given by 
f s with v being the outward pointing normal for all surfaces. Then since ip 
is harmonic we know that 

/ v(ip) + [ + [ = (B.19) 
Js Jt. Jn 
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Now we will calculate what we get for iV and for f u u{jp) for different values of v(if}). 



Set v{ip) = c. Solving the formula for N for |n| we get 

\ 3/2 



N =-\{\L v(v)ll3iA ) 



= -ii' n i 3/2 ^' 2 

|0| = (-4t, 3 ' 2 N/c 2 ) 2 ' 3 
= 2 4 "jrAf 2 <V 4 ' 3 . 

With this in hand we can compute N. 

1/2 



= -& c (X , ' WdA ) 

= _I (22/3^1/2^1/3^2/3) J_ c ^ 

y v(il>)dA 



-2- 1 / 3 vr- 1 iV 1 / 3 c 1 / 3 



If we choose c = —3^7- Then 



AT = -2- 1 / 3 vr- 1 iV 1 / 3 c 1 / 3 ^ 

- 2 - i/3 ^ i/3 (^) v X//H 



B.20) 
B.21) 

B.22) 
B.23) 



B.24) 

B.25) 
B.26) 
B.27) 



B.28) 

B.29) 
B.30) 



Thus m — N is equal to ^ (J* n v(i]j)dA + J s iy(ip)dAy Hence, in order to prove 

monotonicity of X, we only need to show that 47rM > — J e v(tp)dA. 
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Conjecture B.2.1. Let M be an asymptotically flat manifold with positive curvature. 
Let M contain negative point mass singularities Pi and black holes Assume that 
the pi can be resolved with a function if on M so that ip = 1 at infinity, v{p) = 
on E, if = on II, v{ip) = — 071 ^ where N is the mass of the Pi, summed 
appropriately, and II is minimal. Then 

I v 1/2 / v 3/2 

m^ M (M) > ^mflfEi) 2 J - I ^ \m R ( Pl )\ 2/3 1 . (B.31) 

This requires only the following conjecture: 

Conjecture B.2.2. Under the conditions of coniecture \B.2l\ let ip be the harmonic 

function on M so that ip — at infinity, v{iji) = — 1 on S, and ip = on EL Taen 

\ -1/2 / , 
,4 . 4 \ / / 3, 



4tt ( / ip A dAj ip(p 6 dAj > J v{i))dA (B.32) 



Removing this from context gives the following conjecture: 

Conjecture B.2.3. Let M be an asymptotically flat manifold with positive scalar 
curvature. Let M have two sets of boundary £ and II. Let if be a function so that 
ip = on II, v{p) = on £, tp = 1 at 00 and Aip = 0. Furthermore assume that 
v(<p) = 512 ^ji/2 071 n. Tnen Zei ^ 6e the harmonic function defined by ip = on EI, 
ip — ai 00. EFe can choose our boundary conditions for ip on S so £/ia£ 

> ^(^f f A dA\ (^J ip(p 3 dA\ >- J v{i))dA (B.33) 

0/ course this would require that u(tp) > on £ and ?/> < on £. 

If we recast this on the original manifold M, we get the following statement: 
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Conjecture B.2.4. Let (M,g) be an asymptotically flat 3-manifold, with minimal 
boundary E and negative mass singularity p of opposite mass. Also assume that 
(M,g) has non-negative scalar curvature. Assume that p can be resolved by a har- 
monic conformal factor ip so that (M,g) is asymptotically flat, E is still minimal and 
IT, the resolution of p, is minimal as well. We can assume that ip — 1 at infinity, 
<f — only on U, v(ip) = 1 on II, and v((p) = on S. 

Then there are boundary conditions, X , on E so that if ip is the solution to 
A(ip/ip) = 1, with boundary values at oo and p, and X on E, then 

What these X boundary conditions might be we have been unable to ascertain. 
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